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INTRODUCTION 

Calculations  have  been  reported  from  the  Naval  Postgraduate  School  for  the  results  of 
adiabatic  internal  explosions  in  systems  containing  a  variety  of  C-H-N-0  fuels  in  air  with  or 
without  the  addition  of  the  active  metals  magnesium  or  aluminum.1"7  The  present  study 
reflects  the  extension  of  the  studies  to  titanium  as  the  active  metal. 

As  before,  metal  and  fuel  are  considered  introduced  into  a  constant  volume  of  1  m3, 
initially  at  25°C.  Products  are  assumed  distributed  evenly  throughout  the  volume.  The 
process  is  treated  as  adiabatic,  and  the  ideal-gas  approximation  is  assumed  to  hold. 

The  fuels  considered  are  listed  in  Table  1,  along  with  a  few  pertinent  properties. 
Oxygen  balance  is  computed,  on  a  mass-percent  basis,  as  the  excess  (deficiency,  when 
negative)  of  oxygen  in  the  fuel  relative  to  the  production  of  water,  carbon  monoxide,  and 
elemental  nitrogen. 

In  all  the  following,  the  symbol  C  will  represent  the  charge,  or  mass  of  fuel,  per  cubic 
meter  of  air  and  M ,  the  mass  of  metal  (titanium)  per  cubic  meter.  Charge-to-metal  (C/M) 
ratios  and  total  concentrations  (C  +  M,  in  kilograms  per  cubic  meter)  were  varied  over  the 
range  0.1  to  10  for  each. 

In  several  of  the  tables,  a  code  is  used  to  indicate  the  condensed  phases  present.  L  = 
liquid  oxide  solution;  C  =  TiC;  N  =  TiN;  G  -  graphite;  A  =  Ti305;  B  =  Ti02. 


1  Naval  Weapons  Center.  Peak  Overpressures  for  Internal  Blast,  by  G.  F.  Kinney,  R.  G.  S.  Sewell,  and  K.  J. 
Graham.  China  Lake,  Calif.,  NWC,  June  1979.  (NWC  TP 6087,  publication  UNCLASSIFIED.) 

*  Naval  Weapons  Center.  Reactive  Metals  in  Internal  Explosions:  The  Combustion  of  Magnesium  in  Air,  by 
R.  A.  Reinhardt.  China  Lake,  Calif.,  NWC,  February  1978.  (NWC  TM  3429,  publication  UNCLASSIFIED.) 

3  Naval  Weapons  Center.  Adiabatic  Computation  of  Internal  Blast  for  Magnesium-Cased  Charges  in  Air,  by 
R.  A.  Reinhardt.  China  Lake,  Calif.,  NWC,  April  1979.  (NWC  TM  3820,  GIDEP  E202-1481,  publication 
UNCLASSIFIED.) 

*  Naval  Weapons  Center.  Adiabatic  Computation  of  Internal  Blast  for  Aluminum-Cased  Charges  in  Air,  by 
R.  A.  Reinhardt  and  A.  K.  McDonald.  China  Lake,  Calif.,  NWC,  January  1982.  (NWC  TP  6287,  publication 
UNCLASSIFIED.) 

5  Naval  Weapons  Center.  AWorking  Model  for  the  System  Alumina-Magnesia,  by  R.  A.  Reinhardt.  China 
Lake,  Calif.,  NWC,  May  1983.  (NWC  TP  6433,  GIDEP  E392-0754,  publication  UNCLASSIFIED.) 

6  Naval  Weapons  Center.  Computer  Program  for  Internal  Aluminum-Fuel-Air  Explosions,  by  R.  A. 
Reinhardt.  China  Lake,  Calif.,  NWC,  May  1983.  (NWC  TP  6449,  GIDEP  E413-03 19,  publication  UNCLASSIFIED.) 

'  Naval  Postgraduate  School.  Internal  Explosions  of  Reactive  Aluminum  with  a  PBX  in  Air,  by  R.  A. 
Reinhardt.  Monterey,  Calif.,  NPS,  August  1983.  (NPS-61-83-01 1-PR,  publication  UNCLASSIFIED.) 
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TABLE  1.  Properties  of  the  Fuels. 


Designation 

Name  and  chemical  formula 

AUf ,  298- 
kJ/mole 

Oxygen  balance, 
to  CO  +  H2O,  % 

PETN 

Pentaerythntol  tetranitrate, 
C5H8N4012 

-489  8 

+  15 

NC 

Nitrocellulose,  13.3%  N; 
C6H7N25O10 

-8129 

+  3 

HMX 

sym-Cyclotetramethylene- 
tetramtramine,  C4H3N3O8 

+  77  4 

0 

Pentolite 

50%  PETN,  50%  TNT; 
C6  16H6  25N3  4l08  5 

-97  9 

-5 

Comp  B 

65%  RDX,  35%  TNT; 
Cl  96H2.53N2  22°2  68 

+  11  47 

-9 

TNT 

2,4,6-Tnnitrotoluene, 
C7H5N306 

+  30  7 

-25 

N2H4 

Hydrazine 

-165  0 

-100 

C2H40 

Ethylene  oxide 

-47  29 

-109 

Carbon 

Graphite 

0 

-133 

C6H14 

Hexane 

+  10031 

-242 

Internal  energies  of  formation  are  per  mole  of  formula  indicated. 
Data  were  taken  or  computed  from  Ref.  1 


BASIS  OF  CALCULATIONS 

As  is  pointed  out  in  Ref.  4,  the  adiabatic  restriction  requires  that  a  temperature  be 
found  such  that  the  sums  of  the  internal  energies  of  the  equilibrium  mixture  of  products  at 
that  temperature  must  equal  the  initial  internal  energy  of  formation  of  the  chosen  fuel.  The 
pressure  is  then  computed  from  the  ideal  gas  law,  taking  into  account  the  total  number  of 
moles  of  gas  present  at  equilibrium  at  the  temperature.  The  overpressure  is  found  by 
subtracting  1  bar. 

As  before,4  thermochemical  data  are  represented  by  the  five-parameter  expression 

U  =  Bt  +  B2T  +  B3T2  +  B4lnT  +  B5/T 

where  T  is  the  Kelvin  temperature.   In  Table  2  are  listed  the  internal  energy  parameters  for 
the  25  gaseous  and  15  condensed-phase  species  considered.   Those  for  titanium  species  were 
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TABLE  2.  Internal  Energies  of  Combustion  Products. 

Expressed  as  a  function  of  t  =  T/1000  (in  kilokelvins); 

values  given  in  joules/mole:  U(T)  =  Bi  +  821  +  83x2  +  B^lnx  +  Bs/t 


Substance 

B1 

B2 

B3 

B4 

B5 

Ti 

310931 

11660 

75 

1364 

664 

TiO 

264015 

30151 

272 

-2961 

1040 

T1O2 

-55866 

60797 

-355 

-22628 

-4116 

Ti  + 

1029558 

29682 

-252 

-  49343 

-43618 

Ar 

-3718 

12473 

0 

0 

0 

CO 

-  841 1 5 

31136 

-25 

-6218 

24 

CO  2 

-361013 

56789 

67 

-9201 

1764 

H 

213029 

12473 

0 

0 

0 

OH 

166505 

34786 

-88 

-  20093 

-7842 

H2 

136535 

32890 

313 

-21058 

-9413 

H20 

50545 

63224 

-395 

-45949 

-  16417 

NO 

105501 

30619 

10 

-4353 

631 

N2 

38785 

31651 

-61 

-7831 

-657 

0 

214445 

9604 

298 

4428 

2050 

02 

32105 

29476 

662 

-6265 

-2102 

TiC(l) 

-617142 

62760 

0 

0 

0 

TiC(s) 

-762144 

27917 

6225 

14807 

3981 

Ti(D 

-465953 

35564 

0 

0 

0 

Ti(s) 

-695970 

-10976 

8134 

35009 

11727 

C(s) 

-40329 

22904 

320 

10 

0 

CN- 

102390 

31551 

-69 

29 

0 

CN 

691554 

38828 

740 

-37192 

-  24658 

C2H 

758476 

66613 

239 

-43882 

-  20692 

C2N 

538335 

54371 

-  15 

-1052 

2404 

HCN 

284944 

62760 

-302 

-26216 

-5947 

HNCO 

36816 

83592 

-474 

-28281 

-4750 

HCO 

90276 

56166 

-353 

-  19233 

-2345 

CH20 

54534 

85563 

-609 

-32362 

-  1612 

C2H2 

474890 

92019 

-27 

-42593 

-10539 

C3 

906056 

48907 

282 

-16185 

-3220 

TiO  0) 

-1040198 

66944 

0 

0 

0 

TiO  (s) 

-1011246 

56480 

4163 

0 

0 

Ti203(l) 

-1942190 

156900 

0 

0 

0 

Ti203(s) 

-2041671 

145104 

2725 

-30 

4248 

T13O5  <l) 

-2875209 

234304 

0 

0 

0 

Ti305  (s) 

-2963674 

174570 

16844 

266 

123 

Ti02d) 

-  1395260 

87864 

0 

0 

0 

Ti02  (s) 

-1436271 

62766 

5687 

159 

1079 

TiN(l) 

-770564 

62760 

0 

0 

0 

TiN(s) 

-847427 

38404 

5099 

4361 

812 

Substances  are  gaseous  unless  otherwise  specified. 
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calculated  from  data  in  the  JANAF  tables,8  recomputing  the  data  so  as  to  accommodate  the 
universal  choice  of  titanium  vapor  as  the  reference  state.  Parameters  for  the  other  species  are 
taken  from  Ref.  4. 


EQUILIBRIUM  CONSIDERATIONS 


THE  SYSTEM  TITANIUM-OXYGEN 

Unlike  magnesium  and  aluminum,  the  metals  studied  earlier,  titanium  forms  a 
number  of  solid  oxides.  The  compounds  TiO,  Ti203  and  Ti02  are  well  known,  and  their 
formation  accords  with  the  solution  chemistry  of  this  transition  metal.  The  phase  diagram  of 
the  system  Ti-0  is  reasonably  well  established  up  to  the  melting  point.9-10  Maxima  in 
melting  point  appear  to  correspond  to  Ti20,  Ti203,  Ti30s  and  Ti02  TiO  appears  to  melt 
incongruently.  There  is  extensive  solid-solution  formation  and  considerable  temperature- 
dependent  polymorphic  transition.  Oxygen  is  soluble  in  titanium  metal  to  the  extent  of  about 
25  atom-percent,  yielding  solid  solutions  which  show  metallic  conductance.  In  general,  the 
lower  oxides  are  all  nonstoichiometric  and  show  semiconductor  properties. 

The  phase  diagram  gives  no  direct  evidence  regarding  the  liquid  phase,  but  the 
appearance  is  not  inconsistent  with  that  of  a  system  with  no  discontinuities  in  the  liquid 
phase. 

For  these  reasons,  it  has  been  assumed  in  the  present  study  that  the  liquid  phase  in 
the  system  Ti-0  is  continuous;  that  is,  that  there  is  complete  miscibility  from  the  metal  to  the 
composition  Ti02.  This  liquid  is  assumed  to  be  an  ideal  solution  of  the  liquids  Ti,  TiO,  Ti203, 
TisOs  and  Ti02  Ti20  was  not  included  since  no  thermochemical  or  equilibrium  data  appear 
to  be  available  for  it. 

The  complexities  of  the  solid  phases  have  been  ignored  in  the  present  study.  Only  a 
small  number  of  points  were  at  temperatures  below  the  melting-point  curve;  most  of  these 
were  for  oxygen-rich  systems  in  which  solid  Ti02  formed  or  for  carbon-rich  systems  in  which 
TiC  but  no  oxide  at  all  resulted.  Pure  solid  phases  were  assumed  in  those  few  cases  where 
solid  Ti305,  alone  or  with  Ti02,  appeared. 

In  the  computer  program  the  presence  of  "liquid  oxide"  (that  is,  the  liquid  solution  of 
oxygen  in  titanium)  required  in  part  that  the  temperature  be  above  2143  K  (the  melting  point 
of  Ti02  and  the  highest  temperature  on  the  melting-point  curve).  This  criterion  was  adopted 
for  the  sake  of  simplicity,  and  proved  satisfactory  in  all  but  a  handful  of  cases  where  a  liquid 
was  obtained,  but  below  2143  K.  These  points  were  recomputed,  now  waiving  the 
temperature  requirement  for  liquid  oxide.  It  was  found  that  this  change  made  very  little 
difference  (at  most  a  few  tenths  of  a  bar)  in  the  computed  pressure. 


8  National  Bureau  of  Standards.  JANAF    Thermochemical  Tables,  2nd  edition,  by  D.  R.  Stull  and  H. 
Prophet.  Washington,  D.C.,  NBS,  June  197 1 . 

9  R.  C.  DeVries  and  R.  Roy,  Am.  Ceramic  Soc.  Bull.,  Vol.  33  (1954),  pp.  370-72. 

10    P.  G.  Wahlbeck  and  P.  W.  Gilles,  J.Am.  Ceramic  Soc,  Vol.  49  ( 1966),  pp.  180-83. 


NWC  TP  6544 


A  set  of  calculations  was  carried  out  for  the  system  TNT-Ti-air  using  the  quite 
different  (and  less  reasonable)  assumption  that  the  liquid  oxide  phases  were  completely 
immiscible  with  one  another  and  with  the  liquid  metal.  Results  were  appreciably  different,  as 
to  both  adiabatic  temperatures  and  pressures  found.  The  pressures  tended  to  be  higher  (by  up 
to  12%)  even  though  the  temperatures  were  lower.  This  result  was  in  accord  with  the 
reduction  in  fugacities  accompanying  solution  formation. 

It  should  be  remarked  that  it  was  assumed  that  TiC,  TiN  and  graphite,  when  present, 
existed  as  pure  phases. 


EQUILIBRIUM  DATA 

As  in  the  previous  study,4  equilibrium  constants  of  formation  were  fitted  to  the  four- 
parameter  equation: 

log10K  =  At+  A2/(A3  +  T)  +  A4T 

In  Table  3  are  listed  the  equilibrium  constant  parameters  for  all  the  species  considered 
(except  for  the  elements  in  the  reference  states,  for  which  all  the  parameters  are  zero).  Those 
for  the  titanium  species  were  computed  from  the  JANAF  tables,8  with  the  data  below  3591  K 
(the  normal  boiling  point  for  titanium)  recomputed  to  take  into  account  the  choice  of  titanium 
vapor  as  the  reference  state.  (This  choice  was  made  to  avoid  the  discontinuities  that  would 
otherwise  occur  at  the  melting  and  boiling  temperatures  of  titanium.)  Parameters  for  the 
other  species  are  taken  from  Ref  4. 

In  application,  each  Kp  (expressed  in  partial  pressures)  is  first  converted  to  a  Kn  (in 
terms  of  mole  numbers).  Then  for  any  species  the  number  of  moles  is  given  as  a  known 
function  of  Kn  and  the  "master  variables":  X  =  VO~2'>  Y  =  VTT2;  Z  =  VNJ?;  Ti;  and  Ace,  which 
is  the  activity  of  carbon  (standard  state,  graphite).  (Formulas  in  the  last  sentence  refer  to 
number  of  moles.) 


EQUILIBRIUM  CALCULATIONS 

For  argon,  the  mole  number  is  always  0.4036,  the  number  of  moles  in  1  m3  of  air  at 
298  K,  1  bar.  For  each  of  the  five  remaining  elements,  a  material  balance  equation  may  be 
written  representing  conservation  of  the  number  of  moles  of  atoms  of  each  element.  The  mole 
number  of  each  of  the  chemical  species  present  is  a  function  of  those  of  the  elements  in 
standard  states  and  the  activity  of  carbon;  hence,  Five  simultaneous  equations  in  five 
unknowns  (the  five  "master  variables"  referred  to  earlier)  are  obtained. 

In  actuality  the  maximum  number  of  unknowns  to  be  considered  is  four,  since  Y  = 
VH2  can  always  be  found  in  closed  form  as  a  function  of  the  others.  Considering  also  the 
phase  rule  restriction,  the  number  of  variables  to  be  solved  for  is  four  less  the  number  of 
condensed  phases.  The  Newton-Raphson  method,  as  described  in  Ref.  4,  is  then  used  to  solve 
this  set  of  simultaneous  nonlinear  equations. 


NWC  TP  6544 


TABLE  3.  Equilibrium  Constants  of  Formation  of  Products. 

(Base  1 0  logarithm  of  the  formation  constant  of  the 

indicated  substance,  expressed  as  a  function  of  x  =  T/1000 

(in  kilokelvins):  logirjK  =  A^  +  A2/(A3  +  x)  +  A4x) 


Substance 

A1 

A2 

A3 

A4 

TiO 

-2.256 

23.248 

-0018 

-0.133 

Ti02 

-7510 

37  603 

-0  008 

-0  067 

Ti  + 

4082 

-39359 

0079 

0  027 

CO 

4593 

6.108 

0030 

-0.067 

C02 

0.087 

20  642 

0001 

-0024 

H 

3  132 

-12016 

0019 

0015 

OH 

0  769 

-1  969 

-0016 

-0015 

H20 

-3056 

13.305 

0014 

-0  005 

NO 

0710 

-4  800 

0.008 

-0  009 

0 

3497 

-13  439 

0010 

0006 

TiC(l) 

-6352 

29.181 

-0  002 

0019 

TiC(s) 

-8  124 

33929 

-0  003 

0.118 

Ti(l) 

-6  033 

21  383 

-0  060 

-0  006 

Ti(s) 

-7  550 

24  361 

0  000 

0.165 

CN- 

0826 

0  005 

-2.418 

0  047 

CN 

5092 

-22  218 

-0012 

-0016 

C2H 

6  529 

-24  269 

-0012 

-0.043 

C2N 

7060 

-29  343 

-O000 

-0  066 

HCN 

1  607 

-6  807 

-0012 

-0  008 

HNCO 

-25  740 

96491 

2  563 

2  862 

HCO 

2  221 

1  209 

0.219 

-0.060 

CH20 

-2  036 

6  770 

0018 

-0016 

C2H2 

2  655 

-11  332 

-0  016 

0  000 

C3 

10  850 

-41.395 

-0  016 

-0  186 

TiO  <l) 

-8  512 

44816 

-0  050 

0052 

TiO  (s) 

11.560 

50  682 

-0  009 

0  348 

Ti203(l) 

-22.269 

112  701 

-0054 

0  030 

Ti203(s) 

-28610 

126047 

-0  007 

0  607 

Ti305(l) 

-36762 

182  554 

-0039 

0085 

Ti305(s) 

-45  523 

200  007 

-0  003 

0  984 

Ti02(l) 

-13912 

66.399 

-0  042 

0042 

Ti02  (s) 

-17.167 

73  460 

-0  003 

0  290 

TiN  (1) 

-5  271 

20  076 

-0672 

-0  437 

TiN(s) 

-12.371 

41  581 

-0  007 

0  196 

Substances  are  gaseous  unless  otherwise  specified. 
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An  initial  approximation  scheme  is  available  wherein  it  is  assumed  that  oxygen  is 
taken  up  in  the  order  CO,  TiOx,  H2O,  CO2.  With  insufficient  oxygen  to  proceed  beyond  TiO, 
TiN  is  assumed.  With  excess  carbon,  TiC  or  C  (solid)  is  assumed.  Based  on  these 
assumptions,  initial  values  of  all  the  master  variables  can  be  found  and  the  computation  may 
be  started.  After  the  first  computation  in  a  series  it  is  a  matter  of  operator  choice  whether  to 
repeat  the  initial  approximation  or  to  use  the  values  of  the  master  variables  from  the  last  run. 

In  any  case,  after  convergence  of  the  Newton's  method  of  calculation,  tests  are 
performed  for  the  presence  or  absence  of  each  condensed  phase.  If  a  change  from  those 
assumed  has  occurred  the  computation  is  rerun.  This  procedure  continues  until  all  the  tests 
are  satisfied  and  no  change  in  condensed  phases  is  predicted.  Further  detail  on  the 
computational  method  with  a  copy  of  the  computer  program  is  given  in  Appendixes  A  and  B  to 
this  report. 


RESULTS 

For  all  the  systems  studied  the  following  were  found:  adiabatic  temperature, 
overpressure,  and  yield  of  each  product  in  terms  of  moles  and  also  as  partial  pressures  for  the 
gaseous  products.  Total  concentration  (metal  plus  fuel)  was  varied  regularly  over  the  range 
0.1  to  10  kg/m3  and  the  fuel-to-metal  ratio  likewise  was  varied  over  the  range  0.1  to  10. 

In  Table  4  the  following  data  are  given  for  titanium  plus  air  (in  the  absence  of  fuel): 
overpressures  (bars),  adiabatic  temperatures  (K),  and  product  yields,  expressed  as  mole- 
percent  for  the  gases  and  total  mole  number  for  condensed  phases.  The  composition  of  the 
liquid  oxide  phase  is  given  as  x  in  the  empirical  formula  TiOx  (x  ranging  from  0  to  2). 


TABLE  4.  Combustion  of  Pure  Titanium  in  Air: 
Overpressure,  Adiabatic  Temperature,  and  Product  Yield. 


Property/ 

Concentration 

kg/m3 

system 

0.1 

0.2 

04 

1  0 

20 

4.0 

100 

Overpressure,  bars 

7.3 

9.7 

11.7 

14.3 

18.4 

26.3 

49.8 

Temperature,  K 

2602 

3423 

4024 

4338 

4393 

4542 

4869 

Mole-%: 

Ar 

1.05 

1.08 

1  07 

0.95 

0.76 

0.56 

0.32 

Ti       • 

0.01 

7.22 

30  08 

50.59 

72.46 

TiO 

0.12 

486 

17.42 

9.85 

5  22 

2.15 

T1O2 

0.01 

0.08 

NO 

2.47 

5.05 

409 

0.05 

0.01 

N2 

80.84 

81.46 

81.12 

74.26 

5959 

43.63 

25.06 

0 

0.32 

3  64 

6.54 

0  10 

0.01 

0.01 

02 

15.31 

864 

223 

Liquid  oxide: 

Moles 

2.04 

2.77 

455 

10  40 

20.55 

43.2 

115.0 

x  in  TiOx 

1.99 

1.79 

1.40 

0.91 

0.57 

0.31 

0.12 
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Table  5  gives  overpressures  in  bars,  with  the  temperature  in  parentheses,  for  the 
titanium-air-fuel  mixtures.  Each  portion  of  the  table  is  devoted  to  a  particular  fuel,  and  the 
fuels  are  arranged  in  decreasing  order  of  oxygen  balance,  starting  with  the  most  oxygen-rich, 
PETN,  and  ending  with  the  most  oxygen-deficient,  hexane.  In  each  case,  C  +  M  represents 
the  total  concentration  in  kilogram/meter3,  and  C/M  represents  the  mass  ratio  of  fuel  to 
metal.  For  each  table  entry  a  code  indicates  the  condensed  phases  present;  this  code, 
described  earlier,  is  also  explained  as  a  footnote  to  the  table. 


TABLE  5.  Overpressure  (in  bars),  Adiabatic  Temperature  (in  K)  (in  parentheses), 
and  Condensed  Phases  for  the  Combustion  of  Titanium  With  Fuels  in  Air. 


C  +  M 


C/M 


0.1 


02 


04 


2  0 


a.  Fuel:  PETN. 


40 


100 


0  1 

7  1 
(2496) 

L 

68 

(2403) 
L 

64 
(2249) 

L 

58 
(2004) 

B 

5  1 

(1743) 

B 

44 

(1527) 

B 

38 

(1344) 
B 

02 

96 

(3314) 
L 

95 
(3228) 

L 

93 

(3095) 

L 

8.9 

(2832) 

L 

82 

(2564) 
L 

7.5 
(2286) 

L 

67 

(2041) 

B 

04 

11.9 

(3947) 

L 

12.1 
(3874) 

L 

12  2 

(3739) 

L 

12.1 

(3450) 
L 

11.8 
(3209) 

L 

11.5 

(2995) 

L 

11.0 
(2776) 

L 

1  0 

15.4 

(4363) 

L 

163 
(4335) 

L 

17  5 
(4253) 

L 

19.1 
(4033) 

L 

19.6 
(3788) 

L 

196 

(3554) 

L 

194 

(3349) 

L 

20 

207 

(4438) 

L 

22.6 

(4486) 

L 

256 

(4549) 

29.7 

(4349) 
L 

31  6 
(4107) 

L 

32  3 

(3861) 

L 

324 
(3635) 

L 

40 

31.3 

(4586) 

L 

353 

(4639) 

L 

41  7 

(4743) 

L 

51  4 

(4674) 

L 

558 

(4390) 

L 

57  7 

(4116) 

L 

58  3 

(3865) 

L 

100 

63.2 

(4913) 

L 

74.4 
(4970) 

L 

91  9 
(5098) 

L 

121  0 
(5177) 

L 

132.1 

(4773) 

L 

1364 

(4426) 

L 

138  0 

(4127) 
L 

L  =  liquid  oxide,  B  =  T1O2  (s) 
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TABLE  5.  (Contd.). 


C  +  M 


0.1 


cm 


0.2 


0.4 


1.0 


2  0 


4.0 


10.0 


b. 

Fuel:  Nitrocellulose 

0.1 

7.1 
(2508) 

L 

6.9 

(2427) 
L 

6.6 
(2294) 

L 

6.1 

(2089) 

B 

5.4 

(1863) 

B 

4.9 

(1676) 

B 

4.4 
(1  520) 

B 

0.2 

9.6 

(3322) 

L 

9.5 

(3244) 

L 

9.4 

(3125) 

L 

9.1 

(2903) 

L 

8.7 

(2696) 

L 

8.2 

(2488) 

L 

7.6 
(2285) 

L 

0.4 

11.9 

(3949) 

L 

12.1 
(3880) 

L 

12.3 
(3756) 

L 

12.3 

(3497) 

L 

12  2 
(3289) 

L 

12.0 

(3118) 

L 

11.8 
(2965) 

L 

1.0 

154 
(4347) 

L 

16.4 
(4324) 

L 

17.7 

(4246) 

L 

19.4 
(4045) 

L 

202 
(3829) 

L 

20.4 
(3625) 

L 

20.4 

(3451) 

L 

2.0 

20.8 

(4408) 

L 

22.8 

(4428) 

L 

25.8 

(4467) 
L 

30.3 

(4339) 

L 

32.5 

(4120) 

L 

33.6 

(3902) 
L 

34.0 
(3703) 

L 

4.0 

31  3 

(4550) 

L 

35.4 

(4568) 

L 

41.8 
(4610) 

L 

52.9 
(4649) 

L 

57.8 
(4389) 

L 

602 

(4137) 
L 

61.2 
(3908) 

L 

100 

63.0 

(4870) 

L 

73.9 
(4884) 

L 

91.2 

(4922) 

L 

122.9 
(5040) 

L 

137.9 
(4775) 

L 

143.0 

(4436) 

L 

145.2 
(4150) 

L 

c.  Fuel: 

HMX. 

0.1 

7.0 
(2492) 

L 

6.8 
(2396) 

L 

6.4 

(2235) 

L 

(1977) 
B 

49 

(1706) 

B 

4.3 

(1481) 

B 

36 

(1290) 

B 

02 

9.6 

(3311) 
L 

95 
(3222) 

L 

9.3 

(3084) 

L 

8.8 
(2808) 

L 

81 

(2519) 

L 

7  3 

(2218) 

L 

6.5 

(1956) 

B 

04 

11.9 

(3944) 

L 

12.1 

(3869) 

L 

12.2 
(3731) 

L 

12.1 

(3431) 

L 

11.8 

(3175) 

L 

11.4 
(2943) 

L 

10.7 
(2697) 

L 

1.0 

154 

(4340) 

L 

16.3 
(4313) 

L 

17.6 

(4225) 

L 

19.1 
(3999) 

L 

19.6 
(3744) 

L 

19.5 

(3491) 

L 

19.2 

(3267) 

L 

2  0 

20.8 

(4396) 
L 

22.7 
(4406) 

L 

25.7 

(4429) 

L 

29.9 

(4287) 
L 

31  7 
(4032) 

L 

32.2 
(3765) 

L 

32.0 

(3514) 
L 

40 

31  2 

(4536) 

L 

35.2 
(4540) 

L 

41.4 
(4560) 

L 

52.1 

(4579) 

L 

56.2 
(4285) 

L 

57  5 

(3979) 

L 

57.4 
(3691) 

L 

10.0 

62.8 

(4852) 

L 

73.3 
(4849) 

L 

89.9 
(4857) 

I           L 

119.9 

(4910) 

L 

133.7 
(4639) 

L 

1357 

(4236) 

L 

1350 
(3879) 

L 

L  =  liquid  oxide,  B  =  T1O2  (s). 
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TABLE  5.  (Contd.). 


C  +  M 


0.1 


CIM 


0.2 


0.4 


1.0 


2.0 


4.0 


10.0 


d.  Fuel:  Pentolite. 


0.1 

7.1 

(2510) 

L 

6.9 

(2431) 

L 

6.6 

(2301) 
L 

6.1 

(2102) 

B 

5.5 
(1881) 

B 

5.0 

(1698) 

B 

4.5 

(1545) 

B 

0.2 

9.6 
(3324) 

L 

95 

(3247) 

L 

94 

(3130) 

L 

9.1 

(2914) 

L 

8.7 

(2715) 

L 

83 

(2517) 
L 

7.7 

(2323) 

L 

0.4 

11.9 

(3951) 

L 

12.1 
(3883) 

L 

12.3 
(3762) 

L 

123 

(3505) 

L 

12.2 

(3297) 
L 

12.1 

(3130) 

L 

11.8 

(2983) 

L 

1.0 

15.4 

(4334) 

L 

16.4 
(4313) 

L 

17.6 

(4236) 

L 

19.4 
(4037) 

L 

20.1 
(3822) 

L 

20.3 

(3613) 

L 

20.3 

(3432) 

L 

2.0 

207 
(4385) 

L 

226 
(4384) 

L 

256 

(4389) 
L 

30.3 

(4311) 
L 

32.5 
(4091) 

L 

33.4 
(3863) 

L 

33  6 

(3647) 
L 

4.0 

31.1 
(4523) 

L 

350 

(4514) 

L 

41.1 

(4508) 

L 

520 

(4517) 
L 

57.8 
(4343) 

L 

598 

(4071) 
L 

60.3 

(3811) 

L 

100 

62  3 
(4839) 

L 

72.4 

(4822) 

L 

88.3 
(4800) 

L 

117.3 
(4778) 

L 

1378 

(4710) 

L 

142.1 

(4345) 

L 

142.1 

(4000) 

L 

e.  Fuel: 

Comp  B. 

0.1 

7.1 
(2503) 

L 

69 

(2418) 

L 

65 

(2277) 
L 

60 
(2056) 

B 

53 

(1817) 
B 

4.7 
(1618) 

B 

42 

(1452) 

B 

0.2 

9.6 

(3319) 

L 

95 

(3238) 

L 

9.4 
(3113) 

L 

90 

(2876) 

L 

85 
(2648) 

L 

79 
(2414) 

L 

73 

(2188) 

L 

04 

11.9 

(3948) 

L 

12.1 
(3877) 

L 

12.3 
(3750) 

L 

12.3 

(3477) 

L 

12.1 

(3252) 

L 

11  8 
(3065) 

L 

115 

(2893) 

L 

1.0 

15.5 

(4326) 

L 

16.4 

(4302) 

L 

17  7 

(4220) 

L 

19.3 
(4009) 

L 

20.0 
(3776) 

L 

20  0 

(3546) 
L 

19.9 
(3344) 

L 

20 

20.8 
(4373) 

L 

22.7 

(4362) 

L 

25.6 
(4350) 

L 

30.3 

(4270) 
L 

32  2 
(4029) 

L 

328 

(3771) 

L 

327 
(3520) 

L 

40 

31.1 
(4509) 

L 

35.0 

(4488) 
L 

40.9 

(4460) 

L 

51.4 

(4421) 

L 

57.4 

(4264) 

L 

586 

(3949) 

L 

58.2 
(3632) 

L 

10.0 

62.2 

(4822) 
L 

72.2 

(4791) 

L 

87.5 

(4743) 
L 

114.9 
(4655) 

L 

135.4 
(4567) 

L 

139.0 

(4187) 

L 

135.6 

(3756) 

L 

L  =  liquid  oxide,  B  =  T1O2  (s). 
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TABLE  5.  (Contd.). 


C  +  M 


0.1 


CIM 


0.2 


0.4 


1.0 


2.0 


4.0 


10.0 


f.  Fuel:  TNT. 


0.1 

7.1 

(2524) 

L 

7.0 

(2459) 

L 

6.7 
(2352) 

L 

6.2 

(2147) 

L 

5.9 

(2013) 
B 

5  5 

(1862) 
B 

5.1 

(1736) 

B 

0.2 

96 

(3334) 
L 

9.6 

(3265) 
L 

9.5 

(3163) 

L 

9.3 
(2982) 

L 

9.1 
(2832) 

L 

8.9 
(2695) 

L 

8.6 
(2565) 

L 

04 

12.0 

(3954) 

L 

12  2 
(3891) 

L 

12.4 
(3780) 

L 

12.6 
(3547) 

L 

125 
(3355) 

124 
(3205) 

L 

123 
(3082) 

L 

1  0 

155 
(4298) 

L 

16.4 

(4264) 

L 

17.8 

(4195) 

L 

196 
(3995) 

L 

20  4 

(3776) 

L 

20  5 
(3537) 

L 

20  2 
(3292) 

L 

20 

20  7 

(4331) 
L 

22  5 
(4281) 

L 

25  2 

(4201) 

L 

29.5 

(4037) 
L 

32  3 
(3874) 

L 

33  3 

(3641) 
L 

31  8 

(3252) 

L 

40 

30  8 

(4463) 
L 

34.3 

(4404) 

CL 

39  4 

(4352) 
CL 

48  1 

(4194) 

CL 

53.6 

(3979) 
CL 

56 .1 

(3682) 

CL 

54.8 

(3259) 
CL 

100 

60.7 

(4812) 

CL 

69  2 

(4791) 

CL 

82.7 

(4732) 
CL 

1060 

(4531) 

CL 

120  7 
(4244) 

Cl 

1262 

(3833) 

CL 

122.3 

(3286) 
CL 

g.  Fuel:  Hydrazine. 


0.1 

7  2 
(2510) 

L 

7  1 

(2437) 

L 

6.9 
(2320) 

L 

6.6 

(2143) 

D 

6  1 

(1961) 

B 

5.8 

(1809) 
B 

54 

(1684) 

B 

02 

9  7 
(3303) 

L 

9.8 

(3227) 
L 

99 

(3118) 

L 

10.0 
(2921) 

99 

;2?50) 
L 

9  7 
(2589) 

L 

94 
(2436) 

L 

04 

12  3 

(3899) 

L 

12  7 
(3805) 

U 

13  3 

(3659) 

L 

140 

(3386) 

L 

142 

(314'; 
L 

140 
(2895) 

L 

135 
(2644) 

L 

10 

17,0 

(4227) 

L 

188 

(4146) 
L 

21.2 

(3983) 

L 

23  4 
(3555) 

L 

22  8 

(3005) 

L 

20.8 

(2457) 
L 

183 

(1979) 
A 

2.0 

23  4 
(4252) 

L 

27  0 

(4152) 

L 

320 
(3991) 

L 

390 

(3626) 

L 

39  1 

(3042) 
L 

29  9 

(2047) 
A 

28  2 

(1739) 

AB 

40 

36  2 
(4373) 

L 

43  3 

(4248) 

L 

53.4 

(4104) 

NL 

680 
(3735) 

NL 

70  9 

(3101) 

NL 

55  3 

(2062) 

NL 

43  7 

(1443) 

AB 

100 

756 

(4650) 

L 

93.8 

(4482) 
L 

119.9 

(4374) 
NL 

1579 

(3927) 

NL 

165  0 

(3148) 
NL 

130  3 

(2068) 

NL 

93  1 

(1291) 
AN 

L  =  liquid  oxide,  A  =  T13O5  (s),  B  =  T1O2  (s),  C  =  TiC.  N  =  TiN 
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TABLE  5.  (Contd.). 

C  +  M 

CIM 

0.1 

0.2 

0.4 

1.0 

2.0 

40 

100 

h. 

Fuel:  Ethylene  oxide. 

0  1 

7.3 
(2580) 

L 

7.4 

(2568) 

L 

7.5 

(2551) 
L 

7.6 

(2523) 

L  - 

7.8 

(2502) 

L 

7.9 

(2487) 

L 

7.9 

(2475) 

L 

02 

9.8 
(3354) 

L 

100 

(3314) 

L 

10.2 

(3261) 

L 

106 
(3175) 

L 

10.9 

(3103) 

L 

11  1 

(3034) 

L 

112 
(2963) 

L 

04 

12.3 
(3929) 

L 

128 

(3857) 

L 

13.4 
(3743) 

L 

14.3 

(3502) 

L 

145 
(3246) 

L 

142 
(2975) 

L 

137 

(2711) 

L 

1.0 

166 

(4194) 

L 

182 

(4083) 

L 

20  3 

(3903) 

L 

22  7 

(3551) 

CL 

23  5 

(3220) 

CLN 

22  4 
(2760) 

CN 

189 

(2164) 
CG 

2.0 

22.6 
(4207) 

L 

254 

(4085) 

CL 

293 

(3984) 

CL 

350 

(3676) 

CL 

37  3 

(3252) 

CL 

32  9 

(2578) 
CG 

29  6 

(2133) 
CG 

4.0 

34  1 

(4374) 

CL 

39.6 

(4318) 
CL 

47  9 

(4204) 

CL 

599 

(3820) 

CL 

62  6 

(3268) 

C 

57  0 

(2612) 
CG 

51 .1 

(2115) 

CG 

100 

689 

(4741) 
CL 

83  3 

(4676) 
CL 

105  1 

(4535) 
CL 

136  2 

(4022) 

CL 

1368 

(3259) 

C 

129  8 

(2648) 
CG 

1158 

(2107) 

CG 

i.  Fuel: 

Carbon. 

0  1 

7  4 

(2616) 

L 

75 
(2627) 

L 

7  6 

(2643) 

L 

78 

(2667) 

L 

79 
(2678) 

L 

80 
(2678) 

L 

8  1 

(2671) 

L 

0.2 

98 

(3404) 
L 

99 

(3386) 
L 

100 
(3350) 

L 

10  2 

(3241) 

L 

99 

(3038) 
L 

90 

12657) 

78 
(2260) 

L 

0.4 

12.1 
(3979) 

L 

12  3 

(3929) 

L 

12  6 
(3792) 

L 

12.1 
(3336) 

N 

10.5 

(2821) 
CG 

8.4 

(2305) 
CG 

70 

(1960) 

GN 

10 

152 
(4067) 

L 

15.3 

(3908) 

CL 

15  2 

(3804) 

CL 

13.8 
(3588) 

CG 

11.2 

(3010) 
CG 

80 

(2221) 

CG 

58 

(1709) 

GN 

2.0 

192 

(4168) 

CL 

19  1 

(4126) 

CL 

200 

(4473) 

C 

150 

(3788) 

CG 

11.7 

(3128) 
CG 

78 

(2170) 

CG 

5  1 

(1558) 

GN 

40 

265 

(4414) 
CL 

26  1 

(4375) 

CL 

24.5 

(4716) 

C 

16.4 

(3962) 

CG 

12.1 

(3212) 

CG 

7.7 

(2135) 
CG 

3.5 

(1421) 
AGN 

100 

48.2 

(4792) 

CL 

47.8 

(5133) 
C 

36.4 

(5014) 
C 

189 

(4165) 

CG 

12  3 

(3276) 

CG 

76 

(2110) 

CG 

2.4 

(1365) 
AGN 

L  =  liquid  oxide,  C  =  TiC,  G  =  graphite.  N  =  TiN,  A  =  T13O5  (s) 
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TABLE  5.  (Contd.). 


C  +  M 


cm 


0.1 


0.2 


0.4 


1.0 


2.0 


4.0 


100 


j.  Fuel: 

Hexane. 

0.1 

7.6 
(2650) 

L 

7.8 

(2691) 
L 

82 

(2749) 

L 

8.9 

(2824) 

L 

9.3 

(2850) 
L 

96 
(2833) 

L 

9.6 
(2771) 

L 

02 

10  1 

(3387) 

L 

10  4 
(3373) 

L 

109 

(3336) 

L 

114 
(3127) 

L 

11  0 
(2757) 

L 

10.1 

(2391) 

L 

93 
(2082) 

L 

04 

12  7 

(3902) 

L 

13.4 
(3794) 

L 

144 

(3572) 

L 

14.7 

(3037) 

LN 

12.7 

(2359) 

CN 

108 

(1925) 

GN 

9.1 

(1552) 
GN 

1  0 

17  4 

(4005) 

L 

189 

(3805) 

CLN 

20  8 

(3629) 

CLN 

20.2 

(2826) 
C 

17  7 

(2184) 

CG 

15  3 

(1760) 

GN 

126 

(1377) 

ABG 

20 

23  6 

(4075) 

CL 

26.6 
(3972) 

CL 

30.7 

(3740) 
CL 

30  5 

(2891) 
CG 

27  3 

(2162) 

CG 

23  3 

(1674) 

GN 

192 

(1297) 
BG 

40 

358 

(4293) 

CL 

42.3 
(4177) 

CL 

52.2 

(4042) 

C 

51  3 

(2951) 

CG 

46  5 

(2151) 

CG 

39  4 
(1622) 

GN 

31  9 

(1205) 
BG 

100 

734 

(4623) 

CL 

90.2 

(4481) 
CL 

111.5 

(4276) 

C 

113.7 

(3009) 

CG 

^04.3 

(2146) 

CG 

88  2 

(1586) 

GN 

64.2 

0025) 

BGN 

L  =  liquid  oxide,  A  =  T13O5  (s),  B  =  T13O5  (s),  C  =  TiC,  G  =  graphite.  N  =  TiN 


In  Table  6  is  given  the  set  of  product  yields  for  TNT-titanium-air  at  CIM  -  1.0,  chosen 
as  representative.  Mole-percent  in  the  vapor  is  given  for  each  species  that  reaches  at  least 
0.1%  at  some  point.  Numbers  of  moles  of  condensed  phases  are  also  given.  As  with  Table  4, 
the  composition  of  the  liquid  oxide  phase  is  given  as  x  in  the  empirical  formula  TiOK. 

Product-yield  data  are  available  for  all  the  systems  run  but  have  not  been  included  in 
this  report  due  to  the  bulk  of  data  involved.  They  are  available  from  the  author. 


DISCUSSION 

Figure  1  shows  the  effect  of  total  concentration  (C  +  M)  on  adiabatic  temperature  for 
pure  titanium  in  air  as  well  as  for  two  representative  CIM  ratios  for  the  oxygen-deficient 
hexane  and  the  oxygen-rich  nitrocellulose. 
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TABLE  6.  Combustion  of  TNT-Titanium  in  Air  at  CIM  ■  1.0: 
Overpressure.  Adiabatic  Temperature,  and  Product  Yields. 


Property/ 

C  +  M.  kg/m3 

system 

0.1 

0.2 

0.4 

1.0 

20 

40 

100 

Overpressure,  bars 

62 

9.3 

12  6 

196 

29  5 

48  1 

1060 

Temperature,  K 

2147 

2982 

3547 

3995 

4037 

4194 

4531 

Mole%: 

Ar 

0.99 

0.97 

0.88 

0.65 

0  44 

029 

0.14 

Ti 

0.03 

4  47 

680 

784 

TiO 

0  31 

4  35 

3  22 

2.41 

2  57 

Ti02 

002 

002 

Ti  + 

001 

001 

CO 

001 

1.82 

9  79 

2420 

33  46 

3791 

40.95 

co2 

3.78 

558 

3  62 

062 

H 

0  14 

1  53 

7  02 

8  55 

1001 

12  30 

OH 

008 

1  02 

2  19 

091 

001 

0.01 

H2 

0  09 

0  74 

403 

7  65 

1006 

1222 

H20 

1  31 

1  97 

2  19 

0  87 

0  01 

001 

001 

NO 

0  99 

3  08 

3.14 

0  62 

N2 

77.86 

75  58 

69  80 

55  73 

41  72 

31  12 

21  69 

0 

003 

098 

2.66 

0  89 

o2 

14.94 

8  78 

3  13 

0  08 

CN- 

001 

001 

CN 

0  09 

0  25 

0  46 

C2H 

001 

0  06 

0  20 

C2N 

002 

006 

HCN 

0  34 

0.97 

1  52 

HMCO 

001 

HCO 

001 

0  01 

0  03 

0  06 

C2H2 

0  03 

0  11 

Liquid  oxide: 

Moles 

1  04 

1  70 

2.58 

6.48 

13  88 

22  60 

45.39 

x  in  TiOx 

2  00 

1  94 

1  69 

1  18 

0  72 

0  57 

0  56 

Moles  TiC  (I) 

6  17 

29.86 

The  behavior  of  the  pure  metal  is  quite  unlike  that  previously  seen  for  aluminum 
(where  a  maximum  in  T  versus  C  +  M  occurred  at  approximately  the  AI2O3  stoichiometry; 
see  Ref.  4).  With  titanium  there  is  a  steady  increase  in  temperature  with  concentration,  the 
result  of  the  existence  of  a  number  of  oxides,  the  formation  of  each  being  an  exothermic 
process.  Again,  whereas  aluminum  showed  formation  of  A1N  at  high  concentrations,  TiN 
does  not  form  here  as  a  product  owing  to  the  lesser  high-temperature  stability  of  TiN.    The 
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curves  for  the  explosive  (relatively  oxygen-rich)  fuels  are  very  similar  to  those  for  the  metal 
alone. 

With  the  oxygen-deficient  fuels,  first  TiN  and  then  TiC  and  graphite  appear  as 
products  in  the  carbon-rich  region.  These  tend  to  cause  a  falling  off  of  temperature,  as  is  seen 
in  the  case  of  hexane  in  Figure  1. 

Figure  2  shows  the  strong  effect  of  total  concentration  on  overpressure.  Nitrocellulose 
has  been  chosen  as  representative;  all  the  fuels  (except  carbon)  show  a  similar  pattern.  The 
large  increase  in  number  of  moles  of  gaseous  products  combined  with  a  much  less  profound 
change  in  temperature  (Figure  1)  is  responsible  for  the  monotonic  increase  in  pressure. 

The  effect  of  CIM  on  temperature  is  shown  in  Figure  3.  In  general,  titanium  has  a 
higher  heating  value  than  any  of  the  fuels  used;  hence,  there  is  a  generally  observed  decrease 
in  temperature  as  titanium  is  replaced  by  fuel.  At  high  concentrations  with  oxygen-rich 
fuels,  the  oxygen  supplied  by  the  fuel  causes  the  temperature  to  rise  at  first,  resulting  in  a 
shallow  maximum  in  the  curve.  The  considerable  break  downward  in  the  curves  for  hexane 
again  corresponds  to  the  production  of  TiC  and  graphite. 

Commonly,  overpressure  increases  with  CIM,  as  seen  in  Figure  4,  owing  to  the 
increase  in  number  of  moles  of  gases  (CO  +  H2)  as  titanium  is  replaced  by  fuel.  The 
maximum  seen  with  hexane  reflects  the  sharp  drop  in  temperature,  referred  to  above. 

The  nature  of  the  fuel  has  little  effect  on  the  adiabatic  temperature  for  the  oxygen- 
rich  fuels.  For  the  oxygen-deficient  fuels  there  is  a  general  increase  in  temperature  as  the 
oxygen  balance  increases  (becomes  less  negative),  as  is  shown  in  Figure  5,  presumably  due  to 
relieving  the  oxygen  lack.  Pressure  is  surprisingly  insensitive  to  the  nature  of  the  fuel, 
except  for  the  case  of  carbon  which  shows  higher  temperatures  and  markedly  lower  pressures 
as  a  result  of  a  much  smaller  quantity  of  gaseous  product  than  with  the  hydrogen-containing 
fuels.  Carbon  also  shows  a  very  different  set  of  products  from  those  represented  in  Table  6. 
With  carbon  at  high  concentrations  are  seen  large  amounts  of  CN,  C2N  and  C3  and  even 
several  tenths  mole-percent  of  the  ions  as  a  result  of  the  high  temperatures  experienced, 
combined  with  the  absence  of  hydrogen. 


APPROXIMATIONS 

Athow  has  examined11  the  ideal-gas  approximation  for  computation  on  internal 
explosions  in  the  presence  of  the  reactive  metal  aluminum.  In  those  systems  it  appears  that 
the  only  gas  below  its  critical  temperature  is  aluminum  vapor;  and  for  aluminum7  the 
reduced  pressure  can  be  seen  not  to  exceed  1CH4  and  the  reduced  volume  not  to  be  less  than 
200  so  that  no  appreciable  deviation  from  ideality  should  occur.  It  is  anticipated  that  the 
same  conclusions  can  be  made  regarding  the  titanium  systems,  although  no  estimates  of  the 
critical  properties  have  been  made  and  thus  no  estimates  of  errors  attempted. 


11  L.  K.  Athow,  "Real  Gas  Considerations  for  Determining  Physical  and  Thermodynamic  Properties  of 
Gases  Involved  in  the  Prediction  of  the  Effects  of  Internal  Explosions."  M.  S.  Thesis,  Naval  Postgraduate  School, 
Monterey,  Calif.,  June  1982. 
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It  was  pointed  out  by  Smith  that  radiation  effects  should  be  negligible  when 
detonation  takes  place.  When  the  combustion  is  a  deflagration,  however,  the  adiabatic 
assumption  is  expected  to  be  less  reliable.  A  comparison  of  the  experimental  data  on  dust 
explosions  reported  by  the  Bureau  of  Mines13  would  seem  to  show  that  the  computed  results 
from  the  present  report  may  be  high  by  a  factor  of  about  two.  It  has,  however,  been  pointed 
out  by  Baker  et  alM  that  the  Hartmann  apparatus  used  by  the  Bureau  of  Mines  investigators 
seriously  underestimates  the  maximum  overpressures  to  be  expected.  On  this  basis  it  is 
anticipated  that  the  adiabatic  results  will  show  much  less  extreme  negative  departures  from 
realistic  maximum  overpressures. 


12  D.  E.  Smith,  "Attenuation  Effects  of  Thermal  Radiation  on  Internal  Blast  Overpressure."  M.  S.  Thesis, 
Naval  Postgraduate  School,  Monterey,  Calif,  December  1979. 

13  Bureau  of  Mines.  Explosibility  of  Metal  Powders,  by  M.  Jacobson,  A.  R.  Cooper  and  J.  Nagy.  U.S. 
Department  of  Interior,  Washington,  D.C.,  1964.  (RI  6516,  publication  UNCLASSIFIED.  I 

14  W.  E.  Baker,  P.  A.  Cox,  P.  S.  Westine,  J.  J.  Kulesz,  and  R.  A.  Strehlow.  Explosion  Hazards  and 
Evaluation.  Amsterdam,  Elsevier,  1982,  pp.  260-6 1 . 
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Appendix  A 
DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


The  program  itself,  which  is  listed  in  Appendix  B,  is  written  for  the  HP  9845A 
computer.  Various  sections  are  referred  to  by  their  labels. 


MAIN  PROGRAM 

Input  Section 

Fuel  (accessible  from  line  1900).  Enter  formula  and  internal  energy  of  formation. 
Computes  formula  mass;  allows  for  zero  C  or  H  or  no  fuel.  Resets  flags. 

Cone  (accessible  from  line  1870).  Resets  counters.  Enter  C  +  M  (total  concentration 
in  kg/m3)  and  then  CIM,  as  asked  for. 

Temp  (for  entering  total  temperature  manually).  May  be  accessed  by  use  of  special 
function  key  k4.  After  at  least  two  trials,  Temp  may  be  bypassed  and  interpolation  used  to 
find  the  new  temperature.  Pressing  k5,  so  that  At$  =  "Y",  after  two  iterations  allows  for 
automatic  interpolation.  If  needed,  the  automatic  interpolation  may  be  stopped  by  pressing 
PAUSE,  then  k5  (which  returns  At$  to  "N"),  and  finally  k4  for  manual  setting  of 
temperature. 


Computation  Section 

Calc  calls  up  the  computational  subroutines  Eq,  Tical,  and  Ex. 

Eq  is  a  subroutine  of  the  main  program  to  evaluate  equilibrium  constants  of 
formation  of  each  of  the  40  chemical  species.  Kp  is  first  computed  from  the  stored  parameters 
and  then  converted  to  Kn  (expressed  in  mole  numbers).  Kp  or  Kn  is  defined  as  the  ratio  of  the 
activity  of  the  species  to  the  product  of  the  activities  of  C,  Ti  vapor,  H2,  N2,  and  O2,  each  raised 
to  the  appropriate  power.  For  Ti203  and  T^Os,  VKn  is  computed  so  as  to  avoid  overflow  at 
lower  temperatures.  For  condensed  phases,  the  standard  state  is  the  pure  phase;  for  Kn,  the 
standard  state  for  gases  is  1  mole  (in  1  m3). 

Tical  is  the  master  subroutine  which  carries  out  equilibrium  and  energy  calculations. 
Results  are  displayed  as  "dU"  (net)  and  "T  high"  or  "T  low";  dU  =  0  is  desired  for  convergence. 
A  new  temperature  approximation  is  performed  automatically  by  interpolation,  based  on  the 
previous  T  and  dU  values;  or  else  it  is  entered  manually  with  k4.  A  more  detailed  description 
of  Tical  will  be  found  below. 
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Output 

At  Ex,  results  of  the  calculations  are  printed:  temperature,  overpressure,  mole 
numbers  of  all  products  and  the  diagnostic  features  described  under  DIAGNOSTICS,  below. 
Options  are  allowed  at  this  point  for  changing  CIM,  C  +  M  or  for  a  new  fuel. 


SUBROUTINE  TICAL 

Tical  is  the  major  computational  subroutine  whose  task  is  to  find  the  numbers  of 
moles  of  the  products  present  at  equilibrium  at  the  selected  temperature.  The  conditions  to  be 
satisfied  (other  than  for  the  trivial  case  of  argon)  are  the  atom  balances  in  C,  H,  N,  0,  and  Ti 
and  the  establishment  of  chemical  equilibrium  between  each  compound  and  its  component 
elements  in  their  reference  states. 

The  master  variables  (all  in  mole  numbers)  are  X  =  vO~2,  Y  =  VH2,  Z  =  VN2, 
Ace  =  activity  of  carbon  (standard  state  =  graphite),  and  Ti  =  Ti  metal  vapor.  Of  these,  Y  is 
always  computed  in  closed  form;  from  one  to  four  of  the  remaining  are  found  as  unknown 
parameters  in  the  subroutine  Newt  by  using  the  Newton-Raphson  method.  The  actual 
number  of  unknowns  is  equal  to  four,  reduced  by  the  number  of  condensed  phases  present. 
Possible  condensed  phases  are:  Ti  (solid),  TiO  (solid),  Ti203  (solid),  Ti305  (solid),  Ti02  (solid), 
liquid  oxide  solution  (designated  Lox  in  the  program),  C  (solid),  TiC  (solid  or  liquid),  TiN 
(solid  or  liquid).  The  presence  of  each  condensed  phase  is  indicated  by  a  flag,  to  be  set  as 
described  under  Flags. 

Based  on  the  values  of  the  master  variables  and  the  equilibrium  constant  of 
formation,  the  mole  number  of  each  species  is  computed.  Then  the  material  balance  in  the 
elements  O,  N,  C,  and  Ti  is  written  in  terms  of  these  mole  numbers.  When  liquid  oxide  is 
present,  the  stoichiometric  condition  is  that  the  sum  of  the  activities  (that  is,  mole-fractions) 
of  the  components  of  the  solution  should  add  to  unity;  this  condition  replaces  the  atom  balance 
in  titanium  for  this  case.  There  thus  results  a  set  of  up  to  four  simultaneous  nonlinear 
equations;  this  set  is  the  basis  of  the  Newton-Raphson  scheme  to  find  the  unknown 
parameters. 

At  the  conclusion  of  an  iteration,  the  newly  generated  values  of  the  master  variables 
are  used  to  repeat  the  calculations.  In  favorable  situations,  each  iteration  results  in 
improvement  (although  temporary  divergence  sometimes  occurs).  Iteration  is  repeated  until 
the  stoichiometric  errors  are  less  than  one  part  in  ten  thousand. 

Newt  solves  the  set  of  simultaneous  nonlinear  equations  needed  to  find  the  master 
variables.  The  Newton-Raphson  method  is  used.  In  this  subroutine  it  is  necessary  to  find  a 
number  of  derivatives  numerically  by  observing  the  effect  of  a  fractional  change  in  each 
variable.  This  fractional  change  is  set  initially  at  0.1.  It  is  found  that  divergences  which 
would  otherwise  occur  when  the  errors  are  large  can  be  averted  by  altering  this  fractional 
change  to  0.5.  This  is  accomplished  by  the  use  of  special  function  key  k8,  when  the  computer 
is  in  a  pause  mode  or  is  waiting  for  input.  Depressing  k8  a  second  time  will  change  the 
fraction  back  to  0.1.  The  fractional  errors  in  the  stoichiometric  conditions  may  be  observed  by 
the  use  of  kl  (TRACE  VARIABLES  Yn(*)/EXECUTE)  after  which  will  be  displayed  Y(l), 
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Y(2),  etc.,  which  are  these  fractional  errors.  The  display  of  these  will  remain  on  the  screen  if 
PRT  ALL  is  locked  down. 

Approx:  To  begin  the  calculation  an  initial  approximation  of  the  master  variables  is 
needed.  In  this  approximation  an  arbitrary  hierarchy  of  oxygen  and  nitrogen  uptake  is 
assumed.  Oxygen  is  assumed  to  be  taken  up  in  the  order  CO,  titanium  oxides,  H2O,  CO2.  If 
there  is  insufficient  oxygen  to  form  TiO,  then  TiN  is  assumed  to  be  present.  With  excess  C, 
C  (solid)  or  TiC  are  assumed. 

It  has  been  found  that  Approx  gives  quite  satisfactory  values  of  the  master  variables 
at  temperatures  below  about  3000  K;  at  higher  temperatures  the  initial  approximation  may 
be  in  considerable  error.  After  each  entry  into  Temp  the  question  is  asked,  "Do  you  want 
'Approx'?".  The  default  condition  (obtained  by  pressing  CONT)  is  "Y"  for  a  new  concentration 

(C/M  or  C  +  M)  and  "N"  otherwise. 

Flags  sets  the  flags  for  condensed  phases  and  gives  the  values  of  the  temporary 
variables  used  in  Newt.  The  criteria  for  the  presence  of  a  condensed  phase  are  that  the 
quantity  of  the  phase,  if  previously  computed,  be  nonnegative  and  that  the  formation  constant 
be  satisfied  or  exceeded.  In  the  case  of  liquid  oxide,  the  formation-constant  requirement  is 
replaced  by  the  requirement  that  the  temperature  be  above  the  melting-point  curve  (T  = 
2143  is  used  for  simplicity;  see  text  for  discussion  of  this  point)  and  that  the  sum  of  the 
activities  of  the  components  of  the  solution  be  no  less  than  unity. 

Two  important  indexes  are  set  by  Flags.  Iiflag  gives  the  number  of  metal-plus-oxide 
phases  (no  greater  than  two);  for  liquid  oxide  solution,  Iiflag  =  1  except  that,  when  TiC, 
C  (solid),  and  liquid  oxide  are  all  present,  Iiflag  =  0.  The  index  Ii  gives  the  number  of 
unknowns  to  be  sought  in  Newt.  Finally  Flags  gives  initial  values  of  the  variables  to  be  used 
in  Newt. 

Fx  gives  the  fitting  functions  for  Newt.  There  are  three  branches,  depending  on  the 
value  of  Iiflag.  Certain  of  the  master  variables  are  computed  in  Fx.  After  Diff  is  called,  the 
current  mole  numbers  of  all  species  are  generated.  The  subroutine  returns  to  Newt  the 
variable  Fx,  which  is  the  fractional  error  in  the  stoichiometry  in  whichever  element  Newt  is 
considering  at  the  time. 

Diff  first  computes  those  master  variables  which  were  not  found  in  Fx.  Then  Spec 
and  Oxides  are  called  to  compute  the  mole  numbers  of  all  species,  and  finally  the  errors  in 
stoichiometry  for  all  elements  (and  the  sum  of  activities  when  liquid  oxide  is  present)  are 
computed. 

Spec  computes  the  mole  numbers  of  all  gaseous  species,  given  the  current  values  of  X, 
Z,  Ace,  and  Ti.  Y  =  VH2  is  computed  in  closed  form  in  Spec;  it  is  needed  for  computations  on 
hydrogen-containing  species.  For  each  species  the  appropriate  formation  equilibrium 
constant  is  used. 

Oxides  computes  the  mole  numbers  of  condensed  metal  and  oxide  phases.  For  two 
such  phases  (Iiflag  =  2)  the  atom  balance  conditions  for  titanium  and  oxygen  are  used;  for 
Iiflag  =  1,  only  the  titanium  balance  is  used;  for  Iiflag  =  0  the  oxygen  balance  is  used. 
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Nwl  serves  to  route  the  computation  to  the  proper  portions  of  Tical.  On  first  entry, 
the  flags  will  have  been  set  and  control  is  sent  to  Newt;  except,  if  Ii  =  0,  provision  exists  for  a 
closed-form  solution.  In  either  case,  Flags  is  called  again.  If  any  flag  has  changed,  the 
computation  must  be  repeated  for  the  new  set  of  condensed  phases;  otherwise  preparation  is 
made  for  exit  from  Tical  at  the  line  Energy. 

Energy  computes  from  the  stored  parameters  the  molar  internal  energy  of  each 
component  and  then  the  total  change  in  internal  energy  (dU)  from  the  starting  materials, 
taking  into  account  the  mole  numbers  of  each  species  present.  If  the  system  is  at  the  melting 
point  of  any  of  the  condensed  phases  present,  the  amounts  of  solid  and  liquid  are  computed 
from  the  energy  balance,  using  the  known  energy  of  fusion.  Control  is  then  returned  to  the 
main  program. 


DIAGNOSTICS 

Suml  is  called  at  the  end  of  each  run  to  report  the  following  items:  (1)  a  check  on  the 
material  balance  of  each  element;  (2)  the  activity  of  carbon  (Ace)  and  the  sum  of  the  activities 
of  the  components  of  the  liquid  oxide  solution;  (3)  a  comparison  of  the  computed  amounts  with 
the  equilibrium  constants  of  formation  for  each  condensed  phase;  and  (4)  the  activities  (mole- 
fractions)  of  each  component  in  the  liquid  oxide  solution. 

Sum  is  called  whenever  special  function  key  kO  has  been  depressed  once.  At  each 
emergence  from  Newt  the  mole  numbers  of  all  species  are  presented  in  condensed  form;  this  is 
followed  by  Su  m  1 ,  the  output  of  which  has  been  just  described.  These  checks  are  of  particular 
interest  in  troubleshooting.  Depressing  kO  a  second  time  will  cancel  calling  up  Sum 

Depressing  kl  will  cause  execution  of  TRACE  VARIABLES  Yn(*).  Then  during  each 
iteration  in  Newt  the  relative  error  functions  used  to  test  convergence  will  be  displayed. 
Each  Yn(*)  must  drop  below  0.0001  for  convergence.  This  feature  may  be  turned  off  by 
executing  NORMAL. 
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Appendix  B 
PROGRAM  LISTING  FOR  "TIF4" 

10   !  FILE  NAME  'TIF4'  22  Feb  84  :  Ti  t  air  +  fuel;  continuous  liquid  phase 

20   !  Adiabatic  overpressure  calculation  (constant  volume) 

30   ! 

40  OPTION  BASE  1 

50  OVERLAP 

60  ASSIGN  *3  TO  "THTI" 

70  READ  #3,1 

80  DIM  Xn(5),Yn(5),Xt(S),Delx(5),D(5)S)1In(5,5))Nold(S0).Ns(2))Xf(S0) 

90  DIM  N(50),Ak(50,6),K(50),Bu(5D,5),Foi<S0)I8],P(50KNane$(?),U(50) 

100  INTEGER  CflagJinflag.Ticflag 

110  INTEGER  G(50),I,J,K,h,Nw,CuJwct,In,Flag0,Flagl,Flag2,Flag3,Flag4 

120  INTEGER  Ji,IdfJldfuel.Checl()F00,Fi0)F201F3fl.F40)C0)Tin(l)TK0' 

130  MAT  G=CON  !                          Index  for  oases 

140  G(5)=G(13)=G(14)=G(28)=0 

150  FOR  1=1?  TO  23 

160  G(I)=0 

170  KXT  I 

180  FOR  1=35  TO  50 

190  G(I)=0 

200  NEXT  I 

210  READ  *3;Ak(*>, Bud) ,Fo*(*> 

220  LOAD  KEY  "INEK" 

230  DISP  "FOR  ALL  N< I)  PRESS  K0;     FOR  TRACE  VARIABLES  Yn(*> :  Kl" 

240  ! 

250  PAUSE 

260  DEF  FNQd(A,B)C)=(-F+SQR(BA2-4*A»C))/2/A 

270  DEF  FNSti=fi+Tio+Tio2+Tip  !                Sun  qaseous  Ti  atons 

280  DEF  FNSo=Tio+Co+Oh+H2o+No+0+Hnco+Hco+Ch2o+2*(Co2+02+Tio2) 

290  !                                Sub  gaseous  O  atons 

300  DEF  FNOa=Noa-FNSo  !                    Sum  condensed  O 

310  DEF  FNTsu«=Til+2*Ti2+3tTi3+Ti4i             Sun  Ti  oxides 

320  i 

330  C  n=0 

340  Cpn=.l!                              Defaul t  total  nass  is  0,1  kg/cu  n. 


Fractional  change  in  variable  to  compute  derivative 


350  V=l  !                               Default  final  volune  is  1.0  c»  m. 

360  Pv$='V"1                             Default  node  is  constant  volune 

370  Cpn=C  n=.i  !                          Initial  noninal  values 

380  Eps=iE-4  !                           Fractional  error  for  connvergence 

390  Dit=.l  !                             Fractional  change  in  variable  to  c 

400  i                                 Dlt  nav  be  changed  to  or  fron  ,5  by  key  X8 

410  i 

420  Fuel:     INPUT  "Nane  of  conpound",N$ 

430  DISP  "Fornula.  as  atons  per  nole  of  C,  H,  H,  0  in  order"; 

440  DISP  "use  k7  for  TNT"; 

450  INPUT  Ac,Ah,An,Aoa 

460  Met=N(21)=N<22J=Tic=N(i9>=N(20)=Tin=N(43)=N(44)=Gr=N(23)=Cu=0 

470  Cu=Nw=Flag0=Flagi=Flag2=Flag3=Flaq4=Cflag=Tinflag=Ticflag=0 

480  !                                                                                  FIqos  and  condensed  phases  reset  for  new  fuel 

490  IF  Ac=0  THEN  Ac=lE-10 

500  IF  Ah=0  THEN  Ah=iE-10  !                    Provide  for  fuel  with  no  C  or  no  H 

510  DISP  "Internal  energy  of  formation,  J/«ole": 

520  DISP  "use  SHIFT  k7  (=k23)  for  TNT"; 

530  INPUT  Un 

540  Fns=Ac»12.0ii5+Ah*1.008+Antl4.0067+Aoatl5,999i     Fornula  nass,  g/nole 

550  ! 

560  Cone:    OUTPUT  16; "Noninal  C+M  =  ",Cpn 

570  IF  Fns=B  THEN  Fms=1E-10 

580  INPUT  "Total  concentration,  kg/cu  n  of  air,  C+M\Cpn 

590  OUTPUT  16; "Noninal  C/M  =\C  n 

600  INPUT  "Ratio  of  fuel  to  netal,  C/M",C  n 

610  PRINTER  IS  7,1 

620  PRINT  "Total  cone  =  ";Cpn 

630  Mti=Cpn/(i+C  n)!                                                          Total  nass  of  titaniun 

640  Hf=Cpn-Mti  !"                                                               Mass  of  fuel 

650  Cu=Nw=0  !                                                                     Counter  for  iterations  and  criterion  for  'Aoprox' 

660  Anp=PR0UND(Anr3) 

670  Acp=PR0UND(Ac.-3) 

680  Ah  p=P  HOUND  (Ah, -3) 

690  Aop=PROUND(Aoa,-3> 

700  Vp:        Pv»="V 

710  V=i 

720  Cond:              PRINT  USING  Inc;Mti,Mf,N$,Acp,Ahp,Anp,Aop,C  n 

730  Inc:IMAGE  /./,2DZ.3D"  kg  of  Ti",2X,2DZ.3D"  kg  of  "K":  C'^K"  H",K"  N",K"  0",K2X"C/h=  ",3DZ.3D/ 

740  Nti=flti/,0479'                                                                Moles  of  Ti 

750  ! 

i 
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"60  E»i:      Mol5=MAXaDflMMf/F«s.it-iO) 

770  Fns=PROUND(r«s}-3) 

"30  ^o=r^OUND(Mols,5i 

790  Ht=D»0UND(Nti,5) 

S0O  PRINT  «t  ■•  ^ies  9f  Tj» 

310  PRINT  "cormjia  wass  3*  fuel   is":Fr«s  "a/Hole      presen*  are  "'.Mo'.'floles" 

820  PRINT  "VOLUME  =  ';;""  ■:'.!  «" 

830  Conca:  N(6>=  4036  >  A'aon, 

S46  Nc=Ac*riols 

350  Nh=Ah*Hols 

360  Nn=AnMol<;+62^6  '  Includes  noies  of  N  atofts  f-oc  air 

870  Noa=Aoattols+i£.?Sli  '  "    0    "       "       " 

S89  Lf=i  !  Default  value.  TiC  ano  oraohi^e  not  both  present. 

8?0  '_•  D ~ 0  !  Resets  internal  energy 

?00  T i  =U i = 0  '  For  tsflperatwe  interpolation 

?'i0  ftt$=KN"  ;  Default:  narn'ii  temperature  setting 

920  IF  At$0"v"  THEN  Teup 

'30  T=  r590  '•  r "■'  automtic  'efloerature  se*t:na 

940  GOTO  Tenpi 

950 

?60  Teap:     IF  Flacfl  THEN  DISP  "T«=i933  '; 

970  IF  Flacl  THEN  DISP  "Tk=2823  " 

'SO  IF  FIaq2  THEN  DISF  aTw=2ll2  M, 

990  IF  Fiaq3  THEN  DISP  'i^lW  "; 

1000  !F  Flaa4  THEN  DISP  "T*=2143  '; 

iOlG  IF  Ticflaa  THE«J  DISP  NTrt=3290  B: 

1 J 2 Q  IF  linflag  THEN  DISP  hTh=3220  " 

1030  °r eced  1. no  ai»e  «eltina  scints  of  condensed  phases 

1 0 4 C  INPUT  "^pproxi«ate  tewperature  '"',T 

1050  INP'J'  "So  j»oii  wan*  'Approx'W/WJ? 

•.060  !  Default    "Y"  for  -,ew  concentration:   "N"  otherwise 

1070  :  "Approx*  is  Benerally  desirable  below  3000  K,  but  not  above. 

1080  IF  :jfC?.'T$)  =  "Y"  thEN  Nw=S 

1890  Ir  UPC*<T«)="»"  THEN  Na=I 

1180  T$r"     •' 

mo 

1120  Tewpl:! 

1 1 36  Caic:  i>05UB  Ec  '  Climates  eaoiiibriuft  constants  at  T 

ii48  IISr'   "T  =  H-T. 

iiSI  GOSUP  Tkal  :  Subroutine  calculates  loie  nunbe-s  and  '■esiduai  Delta  U. 

il78  D<  =  "  T  n.ah    " 

1180  IF  r=:j'0  THEN  D$="  T  low     " 

1190  DI3T-  M 

1250  :IVFD  0 

;2i9  °PINT  "T  =  ";T;"  K'/DelU  =  ":Dr"  J"  LIN'Z) 

1220  DISP   "T  =  "jT  "  r  SPA(S)   "Deilj  --  •'■[.,   ■  j»; 

12J0  DISP  "« 

1240  STANDARD 

USD  IF  Du=Q  THEN  Ex 

1260  T2=Tt 

1270  U2=UI 

1280  Tl=l 

:290  Ul=Du 

1200  Ir  T2O0  THEN  Tei»2 

1310  IF  At$<>"Y"  THEN  Tewo 

1320  T=3S08-50D*SGN(D«) 

1330  QOTO  TeP)Di 

1340  Tefip2;  IF  ABS!DuM809  THEN  E>: 

13S0  Ir  Att='Y'  "HEN  IntD1  At*  res*t  wi^b  kev  KE 

1360  DISF  "Tr  exit,  press  K2!    nanoai  tewperdtore :  K4'.1 

1270  PIS0  :'  Interpolate:  COM" 

1360  PAUSE 

1390  Intp:  T=T2-U2t(T2-Ti)/(U2-tfi)i  Linear  interpolation 

1400  GOTO  Cak 

tAi  9  1 

1420  Ex:      Nsu«=N*=0 

1430  POP  1=1  TO  44 

14*0  Nt=Nt+N(D!  Total  mies  of  prod^Cs 

1450  Nsu«=Nsun+G(I>W(l!  '  Holes  of  aa; 

1460  NEXT  I 

1470  Conc3;Na=Nsun  '  Moles  of  aos 

1480  HGS5=(ftti+nfH:i*92J*i08&1  Tn»G'  mss  in  arafis 

1490  hasg=f1ass-E9,0itTic-^.9t)iet-.s3  ?«Tii-i47.B*Ti2 

ISOfl  rtasB=«asg-223.7fTi3-79.9»Ti4-61.9itT:n-12,llitGr 

iSifl      !  Kass  or  ]as 
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iS2fl  M  g=Hasg/Ng  !                                                              Avqe  ^olar  miss  of  gases 

1536  M~t=rtass/Nt!                                                                Avqe  nolar  moss  of  all  products 

1541  PRINT  USING  Ihcon;T 

1550  Iwcon:  IMAGE  /*T  =  ",4D.D,"K" 

1560  FIXED  3 

1570  PRINT  "Avq  Molar  MQss(gases)  =  ";M  g;'  g/flol"; 

1580  PRINT  SPA(2)-,"Avg  Molar  nass(total7  =  ":M  f,"  q/hoI"; 

1590  PRINT  SPA(5);Nq:/  Moles  of  GAS  ';Ntj"  doles  TOTAL" ;LIN(i > 

1600  P=8.3143E-5*Ngtt/V           !                                              nPT/v\  VALUE  IN  BARS 

1618  PRINT  LIN(l); 'Overpressure  =  ";P-1;"  &ars";LIN(2) 

1620  STANDARD 

1630  PRINT  •Specie5,;SPA(8):,Moles":SPA(8);,Nole  2";SPA(4);"Part,  pres.': 

1640  PRINT  SPA<4);"LGT(X)";SPA(6);"LGT(P>";LINU> 

1650  !                                                                                 Final  output 

1660  FOR  1=1  TO  34 

1671  IF  G(I)=0  THEN  Conc5 

1680  IF  N(I)(=I  THEN  N(I)=iE-?0 

1690  PRINT  USING  lHCon3;Fo$ d),Ndhl00*Nd)/Nq,Nd)*P/Ng IGTINdJ/NaUGKNdXP/Ng) 

1700  I«con3;  IMAGE  8A,5X,Z.3DE,5X>2D.4D,SX13D.4b,5X)HDZ.3D,5X,HDZ.3D' 

1710  Conc5:NEXT  I 

1720  iKonl:  IMAGE  8A.SX.Z .3DE,5X,Z.3DE, "noles" 

1730  FOR  1=19  TO  23 

1740  IF  Nd)=0  THEN  Conc6 

1750  PRINT  USING  I«conl;Foi(I),N(I) 

1761  Conc6:  NEXT  I 

1770  FOR  1=35  TO  44 

1780  IF  Nd)=0  THEN  Conc7 

1790  PRINT  USING  Inconi;Fo$(I).N(I) 

1800  Conc7:    NEXT  I 

1810  Ceepl:  GOSUB  Sunl 

1820  PRINT  USING  IncjHti ,Mf ,N$,Acp,Ahp,Anp,Aop,C  n 

1830  PRINT  USING  Incf:T,P-i,V 

1840  Incf:  IMAGE  "T  =  '".4D,D\"   Overpressure  =  ',30, 3D."  bars",  "   Volune  =  ",3P.3D,"  cu  n" 

1850  PRINT  LINK) 

1860  BEEP 

1870  LINPUT  "DO  YOU  WISH  TO  RUN  ANOTHER  CONCENTRATION  (Y/N)?\C$ 

1880  PRINT  PAGE 

1890  IF  C$="Y"  THEN  Cone 

1900  INPUT  "Another  fuel?  <Y/N)VCf* 

1910  IF  UPC$(Cf$)='Y"  THEN  Fuel 

1920  OUTPUT  16;  "END  OF  ROUTINE" 

1930  STOP  ! 

i?40  i  tmmmtttmttmmtumttmutumiutmtttutmmmtummt 

1950  Eq:  !     COMPUTES  K.n  FROM  LGT(Kp)  PARAMETERS       tWM**M*tt*mtt»* 

1960  \  mmmtmtmmmmmtmtmmmmtmmtmmttmtmmm 

1970  V  rt=12l88*V/T  '  V/RT  FOR  CONST  VOL  =  V  CU  M. 

l'BO  FUR  1=1  TO  50 

1990  Lgk=Ak(Li)+Ak(L2)/(T+Akd,3))+Ak(I.4'»T 

2000  Lak=Lok+Ak ( I ,6)*LGT ( V  rt)      !                                  Conversion  fron  Kp  to  Kn 

2010  IF  (1536)  AND  (K41)  THEN  Lqk=Lqk/2 

2020  !                                                                                  Use  SOP(K)  for  Ti203  and  T:30S  to  avoid  overflow 

2030  IF  1=43  THEN  Lak=HIN(Lak  .99) 

2040  K(I)=10ALak 

2050  NEXT  I 

2060  Xtic=K(i?>!                                                                For  TiC 

2070  IF  T<3290  THEN  Ktic=K(20) 

2080  Kti=K(2D!                                                                For  condensed  netal 

2190  IF  TU933  THEN  Kti=K(22) 

2100  Ki=K<35)!                                                                     For  TiO 

2110  IF  T(2023  THEN  Ki=K(36) 

2120  K2=K(37)!                                                                  For  Ti203 

2130  IF  T<2112  THEN  K2=K(38) 

2140  K3=K(39>!                                                                     For  Ti305 

2150  IF  T<2047  THEN  K3=K(40) 

2160  K4=K(4D!                                                                    For  Ti02 

2170  IF  T(2143  THEN  K4=K(42) 

2188  Ktin=K(43)!                                                              For  TiN 

2190  IF  T<3220  THEN  Ktm=K(44) 

2208  K043=(K3/K4)A4/Kti/K4 

2210  K032=K2/Kti*(K2/K3*K2t(K2/K3) )'3 

2220  K03i=Kl*(Ki/Kti)t(Ki/Kti)*(Ki/X3)tKi/X3 

2230  !                                These  last  3  used  for  Lox  calcn. 

2240  RETURN  ! 
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2400 


2250  !  tmtmmmummmtmmttmttmtmtmmtmttmmmtmi 

2260  Ticali    PRINT  LIN(2)  !  **»**tt**Wt*m«««*«»**tt*««**«»l***l 

2270  !   WILL  FIND  ALL  HOLE  NUMBERS  AND  THEN  COMPUTE  ENERGY  BALANCE, 

2280  !  mtmtmtmmtnmmmtmtmmMtmtmmttmnmmmtt* 

2290  Cu-0  '                                                                            Cosnts  iiera'ions 

2300  IF  Nw  THEN  Nwl 

2310  Nc:=Nc 

2320  Nnc=Nn 

2330  i  tMtmmutmumttmmtmwmmntxtMummmmmmm 

2340  Approx:  !  (Based  Mostly  on  condensed  phases1   tlt*|*t*$*wt*t*J**l**ll** 

2358  i  wmmmmtmmMtmmmmmmmmtmmmmmtmtm 

2360  FOR  1=1  TO  4E 

2370  N(I)=9 

2360  NEXT  I 

2390  Ar=N(6)=,4036 

400  F:Qa0=Flagl=FlQQ2=Flaq3=F;aq4-Tinflaq=Ticvlaq=Lovfiao=!l 

410  Z=SCR(Nn/2>   !                                                                'Apolies  when  Tinflag=0 

2420  IF  Noa'Nc  THEN  Apox 

2430  Co-Noa  !                                                                             Ng  vide? 

2440  Ticflog=l 

2450  IF  Nti>N;c-Noa  THEN  Adpm 

2*60  Cflaa-1  '                                                                    Graphite  present,  no  Metal, 

2470  Acc=i 

2^99  X=Co/Acc/K(7) 

249(>  Ti=N(i  »=i/Acc/Ktic 

2500  GOTO  Appx 

2510  i 

2528  Appw  Flag  9=1  !                                                                    rte'al  present,  no  trashi^e 

2530  Ti=l/Kti 

2540  Acc=i/Ti/KtiC 

2550  X=Co/X/X(7) 

2560  GOTO  Aop> 

2570  ! 

2580  Aoox:Co-Ncc  !                                                                             Metal  oudes  presen' 

2590  IF  Nti{Noa-Ncc  THEN  ApdI 

2600  Fiaql=l  !                                                                       E*ces;  Ti;   Ti-TiO 

2610  Tinflaa=l 

2620  Exn=Nti-Noa-Ncc-Nn  !                                                     The  excess  of  Ti(c!  aver  N 

2630  IF  Exti>0  THEN  App2 

2640  IF  Ktin/Ki*K2/Ki«2/Kl)i  THEN  Appla 

2650  Z=S8R(-Exti/2)  !                                                          cj  is  in  excess:  TiN  +  T:0 

2u60  Ti=N(i)=i/K  tin/2 

2t?0  ' 

2620  App3;X=i/Ki/Ti 

2690  GOTO  Appx 

2750  i 

2710  ADp2;rlaa0=i              '                                                          Ti(c)  is  in  excess  avtr  N 

2720  Ti=HIN(Exti.i/Kti; 

2,30  Z=i/Ktin/Ti 

2740  GOTO  AppJ 

2?50  i 

2"60  AppliIF  Nti<2t(Noa-Ncc)/3  THEN  Aop4 

2770  X=Kl»Ki/K2/k2        !                                                            TiO      ^203  -<>q;on 

2780  Flaql=Flag2=i 

2790  Ti=l/Kl/X 

2800  IF  Ktin/Kl*K2/Ki*K2/KiU  THEN  Aqpx 

2310  Tinflag=i 

2820  Z=l/Ktln/Ti  !                                                              TiM  2  oxides 

2830  IF  Nti-2*(Noa-Ncc)/3/Nn  THEN  Appx 

2840  If  6*LOG(!f3)+LOG«Ktin);i0*LOG(K2)  THEN  Ano^a 

2850  i 

2860  AoplQ:Z=30R((Nn-Nti+2*iN!ia-Ncci./3)/2)1                                   Ti«  +  Ti203 

2879  Ti=i/Ktin/Z 

2880  X=(l/K2/Ti)*(2/3) 
2890  GOTO  Appx 

2*9  0 
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2910  App4:  IF  Mti< .6*(Noa-Ncc>  THEN  AcpS 

2920  X=(K2/K.3)A4*K2»K2      !  Ti203  -  T:305  rsaion 

2930  Flaa2=Flag3=l 

2940  Ti=f/K2/X/SQR(X> 

2950  IF  6*L0G(K3)+LOG«Krin)<ifl*L0G(K2)  THEN  Aodv 

2960  Tinflag=i 

2970  Z=i/Ktln/Ti    !  TiN  +  2  oxides 

2980  IF  Ntj-3*(NoQ-Ncc)/5>Nn  THEN  Appx 

2990  ' 

3000  App4a:Z=SGR'(Nn-Nti+3*(NoG-Ncc)/5)/2)i  TiN  <■  Ti305 

3010  Ti=i/Ktin/Z 

3020  X=<i/K3/Ti/S9R<Ti})\4 

3030  GOTO  App>: 

3040  ! 

30S0  App5:IF  Nti<(N(w-Ncc>/2  THEN  Aop6 

3060  Floq3=FIaq4=l  '  T =305  -  Ti02  reoion 

3070  X=K3/X4/K4/K4«3 

3080  ! 

3090  App7;Ti=i/K4/X/X 

3100  ! 

3110  Appx:Acc=Co/X/K(7)!  Exit  ffin  "Appro*3 

3120  At3=Ti,3*XAS*K3tK3  '  (Initial  apprtx , ) 

3130  GOSUB  Oxides 

3140  GOSUB  Flags 

3150  GOTO  Nwi 

3160  ! 

3170  App6:Flag4=i  !  Ti02  *  H20 

3180  II7  Noa>lte+2*«ti+Nh72  THEN  Appc 

3190  H2o=Noa-(Nc+2*Nti;   i  H2  -  H20  nixtwe 

3200  H2=Nn/2-H2c 

3210  X=H2o/H2/K(12i 

3220  GOTO  Apn7 

3230  ' 

3240  Appc: If  Noa)2*(lfc+Nti.)+Nh/2  THEN  Aopo 

3250  Co2=Noa-Nc-2*Nti-Nh/2  !  CO  -  C02  wixture 

3260  Co=tfcc-Co2 

32"' 0  X=K(7)tCo2/Co/K(8) 

3280  GOTO  App7 

3290  ' 

3350  Aopo:Co2=Nc  !  Excess  02 

3310  Q2=Noa/2-(Ncc+Nti)-Nh/4 

3320  X=S9R(02> 

3330  Co=K(7)tCo2/X/k<8> 

3340  GOTO  Adp7 

3350  !  tttmmmtmm%tmmttmmmtummmmmtummm%tt% 

3360   '  *  Set-up  for  Newton's  Method  calculation  t 

33?o !  ttmmmmnnimmitmmmmmmmmttmtmmmmmi 

3381  Nwi:  ' 

3390  Nw=Nw+i  !  Controls  entru  into  "Aanrox" 

3400  IF  h  THEN  Nw3 

3418  PRINT  "Closed-fort*  calculation'1  "Newt"  not  used  in  this  case 

3420  GOSUB  Fx2 

3430  GOSUB  Prtflq 

3440  GOSUB  Dspflg 

3450  IF  Check  THEN  GOSUB  Sum 

3460  GOTO  Nw2 

3470  Nw3:  GOSUB  Newt 

3480  Nw2:  GOSUB  Flags 

34^0  IF  (TinfloaOTinO)  OR  (FlagOOFOO)  OR  (FlaoK)FiO)  THEN  Nwi 

3503  IF  <Flag2OF20)  OR  (Flag30F3D)  OR  <Flao40F48>  THEN  Nwi 

3510  IF  (TicflagOTicO)  OR  (CflagOCO)  THEN  Nwi 

3520  IF  (Loxflag=Lox8)  AND  (LflflpLI)  THEN  Ene^av 

3530  I  If  anv  flap  changes,  , tgrations  are  -epeated, 

3540  GOTO  Nwi  ! 
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3550  i       mmummmmmti 

3560  !    *   PRINT  FLAG  STATUS   t 

3570  !       mmmtmmmmm 

3580  PrtflgrIF  Tinfloq  THEN  PRINT  "  TiN"; 

3590  IF  Ticflag  THEN  PRINT  '  TiC"; 

3608  IF  Cflag  AND  (Ac)lE-iO)  THEN  PRINT'  "  Cr": 

3610  IF  Loxflag  OR  Lflag  THEN  PRINT  '  Lox,;,(*;Lf;'),j 

3620  IF  FlagO  THEN  PRINT  '  Met") 

3630  IF  Flogi  THEN  PRINT  "  TiO'j 

3640  IF  Flag2  THEN  PRINT  ■  Ti203"; 

3650  IF  Flag3  THEN  PRINT  '  Ti30S"; 

3660  IF  Flaq4  THEN  PRINT  '  Ti02'; 

3670  PRINT  TAB(i) 

3680  RETURN 

3690  !  mmmmtmttmttmmtmmttmmtmMmmttmtmmt 

3700  i                  DISPLAY  FLAG  STATUS 

37ie  >  mmttmmtmtmtmmmtmtMtmnmttmtmmmmmtmt 

3720  Dspflg:  DISP  TAB(l)   !                       Display  flaq  status 

3738  DISP  "NEWTON'S  CALC,  ';SPA(S); 

3740  IF  Tinflag  THEN  DISP  ■  till"  s 

37S0  IF  Ticflag  THEN  DISP  '  TiC"; 

3760  IF  Cflag  AND  (Ac)lE-iO)  THEN  DISP  '  Gr'; 

3770  IF  Loxflag  OR  Lflag  THEN  DISP  ■  Lox";'(';Lf;*)"; 

3780  IF  FlaaO  THEN  DISP  '  Met"; 

3790  IF  Flagi  THEN  DISP  *  TiO': 

3800  IF  Flag2  THEN  DISP  '  Ti203"; 

3810  IF  Flag3  THEN  DISP  "  Ti30S"; 

3820  IF  Flag4  THEN  DISP  '  TiD2"; 

3830  RETURN 

3840  !  ttmtmmmmmmmtttmmttMmtmmtmmmmmm 

3850  Newt:               !     A  NEWTON'S  METHOD  FOR  UP  TO  FIVE  SIMULTANEOUS  NON-LINEAR  EONS 

3860  i  ttmtmmwmtmmmtmtmtmtmmmmtmmtmmtmmt 

3870  REDIM  Xn(Ii),Yn(Ii)  >Xt(Ii),D€h(Ii),D(Ii>li)  ,In«Ii,Ii) 

3880  GOSUB  Dspflg 

3890  Newt0:Cu=0     !                                                                         Counter  for  Newton's  Method  iterations 

3900  Nwct=Nwct+i!                                                                 Counts  entries  into  "lent 

3910  Newt3;             MAT  Xt=Xn 

3920  PRINTER  IS  7.1 

3930  Cu=Cu+l 

3940  FOR  1=1  TO  Ii 

3958  K=I 

3960  GOSUB  Fx 

3970  Yn(I)=Fx 

3980  FOP  J=i  TO  h 

3990  K=J 

4000  Xt(J)=(i-Dlt)IXn(J)!                                              Dlt=.i  bv  default:  nav  be  reset  r0  .5  by  K8 

4810  GOSUB  Fx 

4020  D(Ij)=-(Fx-Yn(I))/Xn(J)/Dlt!             Partial  of  Ynd)  wrt  Xn(j) 

4030  Xt(J)=Xn(J> 

4040  NEXT  J 

4050  NEXT  I 

4060  IF  DET(D)O0  THEN  Newtl 

4070  ! 

4080  PRINT  'MATRIX  SINGULAR' 

4090  DISP  "MATRIX  SINGULAR' 

4100  PRINTER  IS  0 

4110  PRINT  LIN(i),'D  Matrix;' 

4120  MAT  PRINT  D 

4130  PRINT  LIN(i) 

4140  GOSUB  Sum 

4150  STOP 

4160  ! 

4170  Newtl :MAT  In=lNV(D) 

4180  MAT  Delx=IntYn 

4190  MAT  Xn=Xfi-Delx 

4200  Newt4:F0R  1=1  TO  Ii 

4210  IF  Xn(I))=0  THEN  Newt2 

4220  Xn(I)=(Xn(I)+Delx(I)V2!                                           If  (=  0,  use  1/2  previous  value 

4230  Newt2:    NEXT  I 

4240  FOR  1=1  TO  Ii 

4250  IF  ABS(Yn(I))>Eps  THEN  Neu*3 

4260  GOSUB  Fx 

4270  NEXT  I 

4280  PRINT  Cu;'  ITERATIONS  FOR     ';Ii;'  VARIABLES;  "; 

4290  GOSUB  Prtflg 

4300  PRINTER  IS  7,1 

4310  IF  Check=0  THEN  RETURN 

4320  GOSUB  Sun  !                                                               Diagnostics  given  when  Check=0  (given  bv  kev  K0> 

4330  Retn:  RETURN  ! 
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4340  i  mnmmMttmmttmmttmmmmmmmmwm%mnmm 

4350  Flogs;!  ***********       Sets  flags,  Ii.  initio!  values  for  Newt    #»*****«*»** 

4360  !  ttmtMMtutmmtmtMtiuttmuttmMmMntmnutmmn 

4370  !                                                                                           Sets  flaos  for  condensed  ohases, 

4380  !                                                                                                Formation  constants  bus!  be  exceeded 

42^0  i                                                                                                and  anount  of  phase  not  neoatjve 

4400 

4410  TinO^'inflag  !                                     •                           To  save  old  values 

4420  U=Lflao 

4430  LoxC=LoxfIao 

4440  FOD=FlaqO 

M50  Fl0=Flaal 

4460  F20=Flaa2 

4470  F30=Flag3 

4480  F40=Fl"ag4 

4490  Tic0=Ticflag 

4S00  C0=Cflog 

4510  Cflaa=HlN<i,INT(Acc>)*(Gr>=0> 

4520  IF  GH0  THEN  Gr=0 

4530  IF  Ac-lE-i0  THEN  Cflag=i   !                                              For  fuels  with  no  carbon 

4540  TicflQQ=(TilAcc«tic)I-iE-5)»(Tic>=0)l(Ac)lE-10) 

45'?0  IF  Tk(0  THEN  Tic=0 

4560  Lf=N0T  (Ticflag  AND  Cflaa)1                                          Lf=i  unless  T;C  and  araphite  betft  present 

4570  Tinflac=(Ti*Z*Rtin>i-iE-5)t!un>=0> 

4580  IF  Tm(0  THEN  Tm=0 

4590  FlaoO=<Ti*Kti>i-iE-S?*(lteT>=0) 

4600  IF  He t < 0  THEN  Net=B 

4610  Flagl=(Ti*X«i>i-iE-5)*(Tii>=0) 

4628  IF  TiKO  THEN  Tii=8 

4630  FlBa2=(TitX*SQR(X)$KL'>i-iE-5>*(Ti2>=fl  ' 

4640  IF  Ti2:0  THEN  Ti2=9 

4650  FlaD3=(Ti«MR<Ti)tX»XiS8R(X)IK3>i-iE-S)*ai3>=0? 

4660  IF  Ti3<fl  THEN  Ti3=0 

4670  Flag4=(TilX*XlK4>i-iE-5)*(Ti4>=B) 

4680  IF  Ti4<0  THEN  Ti4=0 

4690  GOSUB  Act 

4700  Loxflag='T>2143)  AND  (Lox>0)  AND  (Ats«n>l-EDs) 

J710  !                                                                                ro-  Lex;  Sim  d*  activities  of  liaoid  oxides  +  netal  )=1 

4720  '                                                                                       and  T  )  neit,  ot.   of  TiD2 

4730  IF  Lox-.O  THEN  Lox=0 

4740  IF  Loxflag  THEN  Flagfl=Flagi=FlagE=Flaa3=Flac4=8 

4750  LfIaq=LoxflQQ  AND  Lf 

4 ''60  IF  Lf lag  THEN  Loxf lac=0 '                                               Lflac  for  Lox  extent  Loxflag  used  when 

477D  !                                                                                              both  TiC  and  graphite  also  present 

4790  Tcon=Nti-FNSti-Tin  '■                                                      Condensed  Ti,  except  for  T;N 

4790  0con=Noa-FNSo     !                                                              Condensed  0 

4800  IF  NOT  (Cf lag  AND  Flaai '  THEN  Flaaa 

4810  Cflaa=<Nce>MlMTcon,Qcon>/3)     !                                  Gr  +  Til    =)  TiC  *  Ti£ 

4820  F:'oQi=l-Cflag 

4630  Flag2=Ticflaa=i 

4340  FlQqa:IF'NOT  (Ticflag  AND  Flag3)  THEN  Flaab 

4850  '  Flag2=Cflag=i  !                                                           TiC  +  Ti3    =  >    C  +  Ti2 

4860  Flag3=(Tcon)i0«Ncc) 

4870  Ticflaq=i-Flag3 

4880  Flaqb:IF  NOT  (Tictlaq  AND  -laa4*  THEN  Flagc 

4890  '  Flaa2=Cflag=i  !                                                             TiC  *  Ti4    =>    C  +  T12 

4900  Ticflag=(210con/'3<Tcon? 

4910  Flag4=l-Ticflag 

4920  Flaqc: IF  <Ticflaq=0>  OR  <FlagO=0)  OR  (Cflaa=0>  OR  <Ac=iE-10)  THEN  Iiset 

4930  '  FlagO-(Nti-FNSti>NccM                                                  Mi*h    TiC.  either  Het  or  G* .  not  Din: 

4940  Cflaa=l-Flaa0 

4950  "  i 

4960  Iiset :Iiflag=Flag8+Flagi+FlaQ2+Flaq3+Flag4+Lflafl  '             hflaa  can  be  0.  i,  or  2, 

4970  '                                                                                                Iifiao  =  i  for  Lflaa  but  0  for  Lexflaa 

4930  IF  hflaa>2  THEN  Flagl 

4990  Ii=4-Iiflag-Tinflaq-TicflQa-Cflaq!                               timber  of  variables  in  *Mevt" 

5900  IF  n=0  THEN  RETURN 

S010  REDIK  Xn(Ii) 

5020  '                                Followmq  sets  values  of  variables  in  "Newt" 

5030  IF  NOT  Tinfiaq  THEN  X»(Ii)=Z 

5040  IF  NOT  Tinflaa  AND  NOT  Cflaq  AND  NOT  Ticflag  THEN  Xn'Ii-i)=Acc 

5050  IF  Tinfiaq  AND  NOT  Cflaq  AND  NOT  Ticfiao  THEN  Xn«Ii)=Acc 

5060  IF  hfiaafe  THEN  XndMI 

5070  IF  (Iiflag=i)  AND  Tinflaa  AND  NOT  FlacB  AND  NOT  Loxflaa  AND  NOT  Lflaa  THEN  Xn(i)=Ti 

5080  IF  (liflag=0)  AND  NOT  (Ctlaa  AND  Tfcflao)  THEN  Xn(2)=Ti 

5090  RETURN 

Siflfi  !                                                                                           The  variables  are  X  «if  no  wore  than  one  oxide  or  ratal), 

5110  !                                                                                               Ti   iif  no  oxide  or  *etai>  and  Z  (if  TiN  is  absent^, 

5l20  i                                                                                           Exception;  When  UN  and  one  solid  oxide,  put  no  petal  or  liauid  oxide 

5130  !                                                                                           present    then  Ti  (and  not  v<  is  first  variable, 

i 
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5140  Flagi;Tcon=Nti-FNSti-Tin-Tic  ■  The  case- of  Ii*laq>2,  which  is  disallowed 

5150  IF  Iiflag=4  THEN  F14 

5160  IF  NOT  FlagO  AND  NOT  Flaqi  THEN  F13a 

S170  IF  NOT  FlagO  AND  NOT  Fla64  THEN  F13b 

5180  IF  NOT  Flaa3  AND  NOT  Flcq4  THEN  F 13c 

5190  IF  Ocon)Tcon  THEN  Fla    !  The  case  of  Iiflaa=5 

5200  Flaq0=Flaai=i 

5210  FlQq2=Fiag3=Flog4=0 

5220  GOTO  Iiset 

5230  Fla;  Flaq0=0 

5240  GOTO  F14a 

5250  F14:  1  For  hflag  =4 

5260  IF  NOT  FlaqO  THEN  F14a 

5270  F13c:  1  Flaas  0,1,2 

5280  Flag2=':Ocon)Tcon' 

5290  Raq>NOT  Flaa2 

5300  GOTO  Iiset 

5310  F13b.  "  Fhas  1,2.3 

5320  FlaaQ=Flaq4=0 

5330  Flag3=(OccnM,6*Tcon) 

5340  FIaql=NGT  Flaq3 

53S0  GOTO  Iiset 

536fi  F14a:   !  '  Iiflaa  =  4.  FlaoO  =  0 

5370  IF  Qeon(=i  5*Tcen  THEN  F13b 

5380  Haal=0 

5390  F13a:  !     '  Flaqs  2,3.4 

5400  FlQg4=(Ocon>5*Tcori/3) 

5410  Flag>NOT  Flaq4 

5420  GOTO  Iiset 

5430  '  ttttmmtimttmmmmmmmmmmtWMmmttmtmmmt 

5440  Fx;!  mtttm      Fitting  functions  for  Newton's  Method  calc,    mttttmtm 

5450  '  tmttmmmMmlmtmMmmmmtmmmmtmttmmtmt 

5460   '  *Fx"  computes  errors  in  stoichiometric  conditions  for 

5470  !  exisitina  "dues  of  variables, 

5480   !  Finds  appropriate  values  of  master  variable* 

5490  1  Calls  or  "Soec"  to  comoute  moles  of  all  species 

5500   '  Xt(j)  is  used,  since  be th  onqinal  value  of  variable 

5510  !  and  value  changed  by  Dlt  are  employed, 

552C  V  NOT  Tin-flag  THEN  Z=Xt(Ii) 

S530  IF  NOT  Tintlaa  AnD  NOT  CHaq  AND  NOT  Ticflaa  THEN  Acc=Xt(Ii-i) 

5540  IF  Tnflaq  AND  NOT  Cflaa  AND  NOT  Ticflaa  THEN  Acc=Xt(Ii! 

5550  ON  Iiflag+i  GOTO  FxO.Fxi..Fx2 

5560  '   ' 

5570  FouTiGOSUP  Diff  ;  Exit  routine  for  "Fx" 

5580    '  Followino  are  stoichiometric  conditions  to  use  in  "New* 

5590  IF  <I=i)  AND  (Iiflao<2)  THEN  F>=Do/Noa 

5600  Ic  (1=1)  AND  (Iifiag=i)  AND  Tinf laa  AND  NOT  Lflog  AND  NOT  flaaO  THEN  -t-^i'^u 

5610  IF  Lonflac  AND  !(>2)  OR  (1=1)  AND'Cfiaa  AND  Ticflaa)  THEN  r^Atsum-i 

5620  IF  '.;=2)  AND  NOT  Iifiaq  AND  NOT  Lflaq  TREN  Fx=Bti/Wi 

5630  IF  WOT  Tinflag  AND  (I=Ii)  THEN  Fx=Deln/Nn 

5640  IF  NOT  Tinflad  AND  NOT  Cflag  AND  NOT  Ticflaq  AND  (I=Ii-i)  THEN  Fx=DeIc/Nc 

5650  IF  Tinflaa  AND  NOT  Cflaa  AND  NOT  TlCflag  AND  (I=Ii)  THEN  Fx=Delc/Nc 

j660  '  The  stoichiometric  teste  correspond  with  ttie  variables, 

5670  RETURN  !  Returns  to   'Newt'  (or  '0  "Nw"  for  Ii=i) 

5680  1 

5690  FxO'    X=Xt(D!  ror  Iiflaa=Q 

5700  IF  NOT  'Cflaq  AND  Ticflag)  THEN  Ti=N(i)=Xti2) 

5710  IF  [flag  AND 'Ticflaq  THE»  Ti=N(i)=i/Ktic 

5720  IF  NOT  Loxflaq  THEN'Fout 

5730  GOSUB  Act 

5740  GOTO  Fout 

5759  >  %mMmmmmmmmm%mttmm%tmtttmmtmt%%ttmm 

5760  Act,'  tmmm    Gives  activities  in  liquid  oxide  tmmmmmmmtt 

5770  1  ttmttMtmtmtmmmmmmmmmtttmmmmnmmtmt 

5780  AtO=Hti*Ti:  Activity  of  metal  in  iinuid  oxide 

5790  '  Followina  a^e  activities  of  oxides  in  hauid  oxide 

5800  Ati=Ki*TitX 

5810  At2=TitTimxtX*K2*K2 

5820  At3=Ti*Ti»TitX»XmX*XtK3*X3*afl0Q<2' 

5830  At4=K4*Ti*X*X 

5840  Atsum=At0+Atl+At2+At3+At4 

58S0  RETURN 

5860  !  mtmmmmummnmmtMmmmtMmmmmmtmtmt 
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5870  F x  1 : '                                                                                     For  Iiflao  =1 

5888  IF  Tinflag  AND  NOT  FlagO  AND  NOT  Ltlaq  THEN  Fxia 

5890  X=X  T ( i  ■ 

S'OO  IF  NOT  Lflag  THEN  Fxib 

5910  !                                                                                For  Lf Iaa=i ,  solves  a  cubic  to  find  Ti  -  f(X') 

5<?26  !                                                                                Must  satisfy  su«  of  activities  =  1 

5930  A3=-i/K3/K3/X/X/X/X/X 

5940  A2=-  ( Kti+KitX+K4*X*X )  *A'i 

5950  Ai=-X*X*X*K2*K2*A3  !                                                  Coeff,  of  cubic  (Unit  coeff.  for  cubic  tern.) 

5969  Ar:(3f  A2-AltAli/3 

5970  Pr='2«Ai*Ai*Ai-?*Ai*A2)/27+A3  !  Reduced  cubic    lacks  quadratric. 
F980  Dc=BrJBr74+Ar*ArtAr/27  •                                            Piscriwinant  of  the  cubic 

5990  IF  Dc<=0  THEN  Fxic 

6000  Ri=-Br/2+SQR(Dc) 

6019  Ri=SGN(Rl)t(SGN(Ri)tRl)"(i/3)  !  Avoids  cube  roots  of  neaative  no. 

6020  R2=-Br/2-SQR(Dc> 

6030  R2=SGN(R2)»(SGN(R2)IR2)Mi/3) 

6040  Xr=Rl+R2  '                                                                  The  only  real  roof  of  the  reduced  ean, 

6050  Ti=N<i)=Xr-Ai/3 

6060  IF  Ti)fl  THEN  Fout 

60?0  fxic      IF  Lflaq  THEN  Ti=N<l)=FN3dvXtXMK2tK2.Kti+Kl*X->-K4MX  -i) 

fc 08**  ;                                                                                           In  case  cubic  had  neaative  •"oot 

6090  Lf=Z 

6100  GOTO  Foot 

6118  Fxlb:     ! 

6120  IF  FlaaO  THEN  Ti=N(i)=i/Kti 

6130  IF  Flaal  THEN  Ti=N(i)=i/Kl/X 

6140  IF  Flag2  THEN  Ti=N(l)=i/K2/X/SQR(X) 

6150  IF  Flag3  THEN  Ti=N(I)=<i/K3/X/X/SQR(X> )*(2/3) 

6168  IF  Flaa4  THEN  Ti=N(i)=i/K4/X/X 

6170  IF  NOT'  (Cflacj  AND  Ticflaa)  THEN  Fo«t 

6130  2=Xt < i > !                                                                     Case  of  TiC  +  araphite  -  Ti2C3 

6190  T|=N(l)=i/Ktic 

b20fl  X=(i/K2/Ti)*(2/3) 

6210  GOTO  fout 

6220  fxia:     Ti=N<l>=Xt(i)     J                                                      Fe:  TiN  ■»■  one  oxide 

6230  IF  Flaal  THEN  X=i/Ki/Ti 

6240  IF  Flad2  THEN  X=(i/K2/Ti>A(2/T 

6250  IF  Flaa3  THEN  X=<i/K3/Ti/S9R(Ti)>\4 

6266  IF  Fia64  THEN  X=S8R(i/K4/Ti> 

6270  GOTO  coi" 

6280  i 

6290  Fx2:      IF  NOT  -FlaaO  AND  Flaal'  THEN  F«2a 

6300  Ti=N(i)=i/Kti  !                                                          ri  -  TiD  region 

6311!  <=1/T:/K1 

6320  GOTO  Fout 

6330  Fx2a:     IF  NOT  (Flaal  ANL-  Flao2)  THEN  Fx2b 

6340  X=Ki/K2*Ki/K2  !                                                           !iO  -  T1203  reoion 

6350  ri=N(i)=i/Ki/X 

6360  GOTO  Fout 

5370  Fx2b:     IF  NOT  (Flag2  AND  Flaq3)  THEN  Fx2c 

63«0  X=<'K2<X3)A4*K2*K2      '                                                   Ti203  -  ^305  reaion 

6390  Ti=*Hi)=i/K2/X/SQR(X! 

640 t  GOTO  Fout 

54)0  fx2c:     IF  NOT  (Flag3  AND  Flaa4)  THEN  Fx2d 

6420  X=K3/K4/X4/X4«K3 

6430  Ti=N(i)=i/K4/X,'X 

6440  GOTO  Fout 

6450  Fx2d:     IF  NOT  (Flag2  AND  Flaq4)  THEN  Fx2e 

6460  X=K2/M*K2/K4 

6470  Ti=Na>=i/X4/X/X 

6480  GOTO  r0ut 

6490  Fx2e:     IF  Cflao  AND  Flaql  AND  NOT  Ticflag  THEN  Alert2 

6500  Ti=l/Xtic 

6510  X=i/Kl/Ti 

6520  GOTO  Fout  ! 
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6330  '  %um%untm%mtntttntmtttt*tmuiw**tt\i%unixmmtnnmt 
6540  *m  Diffi  Calculates  error  function?  for  C  0  N,  Ha  and  ft]  **mt*nf 
isso  >  tmtmtmnmmmimmttttmmtmtmttittitttummtMmu 

6560  Diff :  • 


657  C 

IF  Tinflag  THEN  Z=i/Ti/Ktin 

o580 

IF  Ticfloq  THEN  Acc=i/Ti/Ktic 

6595 

GOSiJB  Spec 

6601' 

Tin=(Nnc-No-2tN2)*Tinflag 

6610 

Tic-(Ncc-Co-Co2)ITicflac!                                          *"ic  and  Ncc  are  corrected 

662C 

1                                                                               :'0*s    a;  anen  _n  "Spec" 

663! 

N(44)=TlnMT<3228 

6640 

N(43^=Tm-N(44) 

6650 

Ni2€)=Tic*<T<3290) 

6660 

N(i9)=Tic-N(20) 

66711 

SOSUB  Oxides  !  GIVES  OXIDES  AND  METALS 

6683 

Dti=Nti-FNSti-Tin-Tic-Met-F'NT£«« 

6698 

I'o=Noi-rNCo-Til-3^Ti2-5*T13-2*T14 

6700 

De'Jn=Nnc-No-2M2-Tin 

6710 

If  Cflac  m  Ticflaa  THEN  Tic=Dti^K 

6?20 

Deic=Hcc-Cc-Co2-Tk' 

4730 

Gr=NCi)=Delc«Cfiao*(Ac>lE-i8) 

6740 

RETURN 

6750 

i 

W,li    l 

ttmmmmmttmimmtmtmtMttmmimtmmuumtmmt 

6770   ' 

mmtmtum     Computes  noles  of  ail  gaseous  species    *«*#*«**?*»*** 

;.760  ' 

«»t**t*xt*«i?«:*t*«!r*i«ittt****t**»*t»*i*tt*tt«*tttt«t*»;««t«$t***»$**t**t* 

6?ef  Spec;!                                                                                    Input;  V    Ti.  Z.   Ac; 

6800 

D2=N<18)=X*2 

6810 

D=«(i7)=s(i7)»X 

S82C 

Tio=N(2)=K(2>«Ti*X 

6830 

"io2=M»  3  "=K  i  j;*T.  *X*X 

3840  N. 

it: 

6850 

N2=N(l6>=Z*2 

bS60 

No=N<15)=K!15)*X*Z 

6870 

rtoot-'=N'42>=K<42^*hci>*X 

icdO  Ci 

:-:     IF  Lfiag=i  THEN  Ace-: 

6890 

I1'  Ac=it-18  Then  Acc=lE-3fl  !                                      p>  fuel*  with  no  carbon 

6900 

Do=N(7)=K(?)*Acc** 

69l8 

Co2=N(8''=k(8)<Acc*X*X 

i920 

Mh=K(?>+K(l0.)*X+Acct<K(2?)*Z+K<30>*Z*XTM3!  *X'  -  Zii?Acc*Ac: 

6930    : 

"'' "* it.-";  t  of  Nonoh"d"efi€ 

S':'fl!) 

Dh=l+Xil2)»X*KiZ2)tAcc*X+K(33i*AcciAcc                       ^efficient  Bf  dihudrooen 

6950 

Y=FNQti(2Wh  Mi.-Nh) 

i960 

H=jHi'9is|((9)tY 

i  a'fi 

Dh=Xi  l  D'-K1 l"i*;V 

svdO 

H2=NQi)=f*Y 

6990 

^2o=*t(i2?=K<i2)*X*¥*if 

7961 

Hcn=*(29'=K<29>*Ao::*Y*Z 

°ui0 

Hncc=Ni30)-K(36)*AcctXtY*2 

7020 

Hcg=N(31)=Ki3i)*Acc*XtY 

7030 

C  rr  2  c  -  N  ( 3  2 .) = K  i  32 .'  *  Ac  c  a  X  *'  V  *  Y 

"940 

C2h2=N  33 )  =X  ( 33 )  *Acc tftc c tt  l ". 

7050 

Tip=N(4)=Cn«=N(24?=S9R(K(4)tK{24)»Ti*Acc?Zi 

7860 

Cn=N(2S»=K(25)lAcciZ 

7070 

C2h=N(26)=K(26)*Acc*Acctf 

"0^0 

C2n=W27)=K<27)*Acc*Acc*Z 

7890 

C3=K(34)*K<34)tAcc*Acc*Acc 

7 1 5  G 

Ncc=Kc-f<cn-Hfico-Hr8-Ch2'!-Cn-Cn«-2*(C2h2+C2h+C2ri   -3*C: 

7110 

Ctr'ecteti  soles  ■'■   C  fTn«; 

'120 

Nnc=Hn-tfcn-Hnco-Cn-C2n-Cn^     '                                      lor-ected  'soles  of  w  a'eits 

7130 

IF  «nc<'>  THE«  Xnc=i£-55 

"140 

RETURN  ! 

coles  e*  w  or  L 


specie; 


38 


NWC  TP  6544 


7iso  !  tmmtmMmttmtmmmtmutmtuttmmmmmtmmmm 

7160  '  ***********  Moles  of  «etals  and  oxides  MM**************************** 

7i'/o  !  tmtmmttmmmmmtmmtmmmmmmMtmtmmmm 

7180  Oxides:  FOR  In=i?  T0  22 

7190  N<>n)=0 

7200  NEXT  In 

?2ifl  FOR  In=35  TO  44 

7220  N(ln)=0 

7230  NEXT  In 

7240  Het=Tii=Ti2=Ti3=Ti4=S 

725(1  '                                                                                Preceding  prevents  carryover  of  old  values 

7260  IF  HOT  Iifiaa  AND  NOT  Loxflag  THEN  RETURN 

7270  Tc8n=Nti-FNSfi-Tin-Tic  !                                           Condensed  Ti  ♦  o> ides 

7280  Dcon=Nca-FNSc     '                                                        Condensed  0 

'290  IF  Lsxflog  OR  Lflag  THEN  Qxlox 

7300  IF  Iiflag=i  THEN  Oxl 

7318  !                                                                                Following  for  two  oxides  or  Metal  +  ine  oxide 

7320  IF  NOT  (Fiaofl  AND  Flag!)  THEN  Oxo 

7330  Tii=Occn 

7340  ltet=lcon-Ti3 

"3b0  GOTO  0x2 

7360  ; 

7370  Oxb;   IF  NOT  (Flagi  AND  ciaq2)  THEN  Oxc 

7380  Ti2^0con-Tcor, 

"390  Tii=Tcon-2tTi2 

V4I0  GOTO  0x2 

7410  ! 

7420  Oxc:      If  NOT  (Flag?  AND  Flaq3!  THEN  Oxd 

?43C'  Ti3=2*Ocon-3tTc8n 

7440  Ti2=CTcon-3*Ti3)/2 

'450  GOTO  0x2 

?4fc3  ! 

7470  Oxd'      IF  NOT  <Fiaq3  AND  F.ca*)  THFn  g>>; 

7480  Ti3<J*Tcon-Ocon 

-490  T.4=Tr0n-3iCTi3 

7500  BOTH  0x2 

7520  Oxe-       IF  NOT  Fiaa2  AND  Flac4  THEN  Alert? 

7330  7i2=2*Tcon-Ocori 

7540  Ti4=Tcon-2»Ti2 

755G  GOTO  0x2 

7560  (hi:      IF  Tin'laq  THEN  Qxiin 

7570  IF  FlacO  THEN  Het=Tcon  i                                            -or  Iiflag  =i.  no  liqind  wide 

7580  IF  Flcci  THEN  Tii=Tcon 

7590  IF  Ficq2  THEN  Ti2=Tc«n/2 

7660  IF  Flag3  THEN  Ti3=Tcon/3 

76i8  IF  Flag4  THEN  Ti4=Tcon 

"620  C0TO  0x2 

7630  ! 

7b40  Qxtin:  IF  Flaafl  THEN  Het=7con-Nnc+No+2*N2    '                      Routine  for  TiN  *  one  ond*.    fhi=  lire  is  mi  bal 

7c;0  Ir  r:aai  THEN  Til=Ocor 

'660  IF  Flag2  THEN  u2=Qcnn/3 

7670  IF  Fiao3  THEN  Ti3=0cxir,/5 

"630  IF  Fiad4  THEN  Ti4=Ocon/2 

7693  GOTO  0x2 

7700 

'710  Qxlox:  S»M9X=Tcon/(H,At2+2tAt3)!                                            Liouid  oxide  oresen* 

7720  IF  Cflag  AND  Ticflag  THEN  B»«ox=Ocon/(i+2*At2+4tAt3+At4i 

7730  !                                                                             For  *hi«  cose  use  0  "o+he''  'ho^  T'  salance 

''740  G0EOF  Ac* 

7750  M.e*=AtO*Simox 

7?60  Tii=AtitSuMX 

"770  Ti2=At2*S«Mox 

7780  Ti3=At3*SuHox 

7790  TiA=At4lSuMox 

'300  Lox=Het*7ii-tTi2+Ti3+Ti4!                                          *c!e?  tf  licud  oxide 

7810  ' 

7820  0x2:      M<22)=Met*(TU933) 

7830  N(2i)=»feT-N(22) 

7840  N<36)=Tiit<T<2023) 

7850  N(35>=Tii~N(3t) 

7860  N(38)=Ti2*(T(2ii2) 

7870  N(37)=Ti2-N(3B) 

"330  H(48l=Ti3t(T<2647! 

7896  N(39'=Ti3-N(4J) 

7V00  N.(42Mi4tfT(2143i 

791C  M(4i?=Ti4-N(42) 

""20  RETURN  ' 
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■■■'70  '  mmmmmtmmtmmmmtmtmtmwmMttmmmmm 

75>40  !  tttmmm         Energy:  Computes  Delta  U.       mtmmmmtmmmw 

?v55  i  ttmnmmmtmmMtmmmtmttummmtmmmmmmtt 

7*60  Ej-ierQvrDu=-US 

7V70  "  FOR  1=1  TO  44 

7980  U<IJ=Bo<I.i)+Bu(1 ,2)*T+Bu(I  .SXTtT^Bu' I  f  5)/T+Bua  .4)*LQG(T) 

7>?H  !                                                                           Molar  internal  enerqu  of  each  specie; 

3000  Du=DuHJ<I)*N(I!!                                                  Softs  total  internal  enerav 

8010  NEXT  I 

8020  IF  Du<=0  THEN  RETURN 

8030  IF  FiaaO  AND  (T=i933)  THEN  EnO 

8C48  Ic  Ticfiaq  AND  • T=32?0 )  THEN  Entic 

8050  IF  Flaqi  AND  (T=2823)  THEN  Eni 

8060  IF  Flaa2  AND  'T=21i2)  THEN  En2 

8070  IF  Flag]  AND  <T=2047»  THEN  En3 

90::0  IF  Flag4  AND  (T=2143)  THEN  En4 

8090  IF  Tinflao  AND  (1=3220)  THEN  En  tin 

81  CO  RETURN  ' 

8110  '                                                                                        In  the  f'oliouina.  2  ohases  of  the  substance  a.--:  present, 

8128  ! 


8130  Entic:  N'2D)=Du/71145  ' 

8140  N<i9)=«AX(i,Tic-N(20))   i 

815G  Di'=Bu*(N(i?)=fl) 

3161  If  N(19)  THEN  RETURN 

8170  N(28)=Tic 

8180  GOTO  Enbeep 

3190  EnO;  N(22)=Du/i8623  !  Soli'  Ti 

8200  N(2i!=MAX(0,Het-N<22>)  '  Liquid  Tj 

9210  Dii=Du*(N«2i)=0) 

8220  IF  M21)  THEN  RETURN 

3230  N<22)=He1 

8240  GOTO  Enbeep 

B2SDEntin:  N(44)=Du/66960  !  Solid  TiH 

8260  N(43'=MAX(0)Tin-N(44)'   !  Licud  TiN 

8270  Du=Du*<N<43)=0) 

3280  IF  N(43)  THEN  RETL'RN 

3290  N'44)=Tin 

&300  GOTO  Enbeep 

8310  Enl:  N(36)=Du/S440S  !  Solid  TiO 

8320  N(35)=HAX(0,Til-N<36))   >  Liouid  TiO 

8338  Du=Du»(N(35J=0> 

334  J  IF  N(35)  THEN  RETURN 

3250  N' 5S)=TAi 

3360  GOTO  Enbeep 

33~0  En2  N(38)=Du/110484  !  Solid  Ti203 

8380  N.(3?)=KAX(0  Ti2-Nl3B))  !  Liaaid  Ti2D3 

3390  Du=Dut(N(37)=fl) 

3401  IF  N(37)  THEN  RLIljR* 

3410  N(38)=TiZ 

8420  GOTO  Enbeep 

8438  En3:  N(4D)=Du/138l05  '  Solid  T>70$ 

8440  N<39>=MAX(0Ji3-H<48>>   ^Liquid  "I30S 

9450  Du=D')t(N(39)=8.) 

8460  IF  Ni39)  THEN  RETURN 

3470  N«40)=Ti3 

8480  GOTO  Enbeep 

8^90  En4-  N«42)=Du/66960   !  Siiid  Ti02 

SS0P  N(4i'=flAX(8,Ti4-*i42.i    1  Liauid  TiQ2 

8bi6  D>j=Dut'.N<4i)=0) 

8520  IF  N(41)  THEN  RETURN 

S':.30  N(42)=Tj4 

8540  ! 

Bb-58  Enbeep  i     8EEP 

3560  DISP  "Select  a  lower  temperature" 

3570  RETURN  ! 


Return 

frofi  "TicaJ 

."  to  nain  p"oa*a« 

In  the 

following 

2  chases  of  +he  subs 

tsnce 

we 

D 

•i 

Solid 

is  ai»en  by 

dU  divided  by  enerav 

e<  f 

iisisr 

Solid 

TiC" 

Liauid 

TiC 
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8580  !  mmmtmmmtmmmmwmnmmmmtumtmmtmtmt 
8590  !  mtmnmmm  Bicqnosucs  tmmmnMttmmmmtnmt 
860 o !'  nnmmmttMmmmMmnmmtmmmmmmtmmmmt 

8610  Alert2:      '  Improper  oxides  present 

6620  PRINT  "Error  in  OXIDES" 

6630  GOSUB  Sim 

3640  STOP 

86S0  GOTO  Tenp 

8660  ' 

8670  ! 

8680  Sun:  PRINTER  IS  7,1!  If  activated  by  Kl,  orints  after  each  iteration 

8690  GOSUE  Spec 

S700  GQSUB  Prtflc 

B?iO       ■     PRINT  LINfD/T  =  M,LINU> 

8720  1=1 

8730  pRIN'  ysiNG  lHSUH;Fo*(I),N<I),Fo*(I+i),N<I+ihFo$<I+2>  NU+2)  .F«*(I+3),N<I+3) 

3740  1=6 

B?5>fl  PRINT  USING  I«stt«tF8$(I)J(I},Fo*(I+i),»«I+i),f<i*(I+25  N(IA2) ,Fo$<I+3) ,«<I+3) 

8760  1=10 

8770  PRINT  USING  1ms'jp-Fo$<I)  ,N(I)  ,Fot!Pi>  N(I+i)  .Fo$(I-2\Na*2>.Fo*<I+3>.N(I+3) 

8780  MS 

B790       PRINT  USING  lHSWfl',Fo$(I),N(I)  JaJd+D.N'.I+i) -.FotKl+2) ,N(I*2),F«*a+3hH(I+3) 

6800        1=24 

8810        PRINT  USING  I«sM'.F»$(I),N<I),F9*(I+i)  N(I+i),F«*(I+2i  N(I+2),Fo$(I+3),N'I+3) 

6820        1=29 

S33C       PRINT  USING  Iksuh;Fo$U)  J(I),Fo$U+l)  N(I*i),F«$U+2)  N(I-,Z)fFi$(I+3)-,N(I+3) 

?O40        '=71 

68S0        PRINT  USIMG  Iw»N}F«*(I),N(I),Fo*a+i)  J(I+i) 

3860  PUNT  USING  Insum'TiC*  Tic, "Met"  .Net  ,'TiOVTii.  "T.2D3"  Ji2 

8870  PRINT  USING  Insun'."Ti30$".Tl3,"Ti02\T.>4  *TiN".Tin  "GrVGr 

8880  Insun:  IrtftGE  3<8A ,fiZ . DD»E, 3a > 

8390  SuhI;  SuM=FNSti+Tin+Tic+Met+FNTsu«-1  Aocears  in  final  output 

8900  PRINT  LINii),"Ti  BALANCE   "Nu  Sum 

8710  3.j«=FN3o+Til+3tTi2+StTi3+2*Ti4 

8'20  PRINT  "0  BALANCE  ".Nea.Sim 

3V30  SuH=Tin+2*N2+No+Hcn+Hnco+Cn+C2n+Cn« 

8940  PRINT  "N  BALANCE  '  Nn,S»fl 

8'50  S'u.h=Nc-Ncc+Co+Co2+Tic+Gr 

B960  PRIM"   '[  balance \ifc,SuH 

897Q  S<JH=«h*Y+2*Dh*W 

8980  PRINT  "H  balance"  .Nh , Sun 

8990  3RINT  USING  "K.K";"Acc=  ".Ace 

7000  PRINT  USING  "K'.K.X/" ■  "Ats»«  =  '.Atsutl 

?8i0  ^RINT  LIN(2>iTABt20) i "Test  for  equilibri!m\LIN(i! 

9B20  PRINT  5PA(  1 ; ,  'K";  rAB'(iO ) ;  "Nuner  ic'jTABGfl )  '■. 

90311  PRINT  "Functional  (<=Nu«e-ic)".LIrJ(i'> 

9940  FLOAT  3 

"OSO  IF  NOT  Loiflag  AND  NOT  Lflaq  THEN  ArO=Ati=At2=At3=A*4=l 

7060  PRINT  "THg}";TAB<iD),fttB.'KtijTAB(33)  Ji 

9070  PRINT  "TiOa;TftB<iO),Ati/Ki!TAB(3(l),TifX 

9080  PRINT  "Ti2O3".TAB(l0))SQR(At2)/K2'.TAB(30))Ti*<iSOR(V 

?CV0  PRINT  "Ti3O5",TAB(l0),SQR(AT3)/K3:TABi30).TitXtX*SOP(Tit>.i 

^100  PRINT  "Ti02b:TABaO'  At4/K4jTAB(3u).Ti*02 

9110  PRINT  ■TiN";TAB<i0U/KUn'.TAB(30>,ti*Z 

9120  PRINT  "Til";TAB(i0),i/!!tic;TAB(30Ui*Arc 

9130  IF  Lo^fiaq  OR  Lflaq  THEN  PRINT  LIN(i) 

9140  IF  Uxflog  OR  Lflac  THEN  PRINT  USING  IissiiH:*a(het>'  AtO  VTiO>"  .Ati  ,"q(Ti2G3)  "  .At? 

MSB  IF  Lcxfiac  OR  Lflaq  THEN  PRINT  USING  lHsur,"a(Ti3Q5J"  At3.  VTt02'\At4 

91.60  !  Actiwitass  in  liquid  sxide 

9170  STANDARD 

9180  PRINTER  IS  7,1 

'190  RETURN  i 
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he 

Ticflaa 

Ktic 

lie 

Ticflaq 

Ktic 

het 

FlaqO 

Kri 

Ket 

FlaqJ 

Kti 

5r 

Cflaq 

■'2oo  <  mtmtwttmmtmmitmxmtmmttmmmmmtMttmtmm 

9210  i  ttmmttltm      Species,  codes,  and  flaas  MMtMttmm**«**mt** 

?22o  \  xttmtmtmmtmmtimmimmmmmmMmmtmmmmt 

9230  '  CODE         rLAC  FORMATION  CONST 

9245  ! 

'250  '  1    h  Ti 

9260  !  2    TiO 

9270  '  3    Ti02 

9230  !  4    Ti- 

9290  i  5  [not  used] 

9390  '  fa    Ar 
"310  '  7    CO 
9320  !  6    C02 
9330  '  Q    H 
9340  '  10     OH 
93S0  '   11    H2 
9360  !  12    H20 

?3?0  '   13  [not  used] 

9380  i  14  [not  used 

9391  '  15  NG 
9400  i  16  N2 
9410  '  17  0 
9420  i  18  02 
7438  i  19  TiC(l) 
9449  !  20  TiC(s) 
«4S8  '21  Ti(l) 
9460  '  22  Tiis.i 
<?470  !  23  C(s) 
°480  !  24  CN- 
9490  i  25  CH 
9500  !  26  C2h 
9'r.lO  i  27    C2N 

9520  I  28  [not  used] 

9'r.30  !  29    HCN 

5548  !  30    UNCO 

9550  !   31    HCO 

9560  !   32    Ch20 

9^0  '  33    C2H2 

9580  !  34    C3 

"590  i  35    TiO(i)  hi         Flea:  Ki 

96  Oil  '  36    TiO(s)  in         Flaqi  Ki 

9610  ■■   37    Ti2Q3(l)  Ti2         Flag?  K2=58R<Xn.' 

9620  !  38    Ti203(s>  Ti2         Fiao2  K2=S&(Kni 

9630  i  39    Ti305U)  Ti3         Flac3  K'3=SQPCKni 

^40  '  4'.'    Ti3G5(5)  ii3         Fiao3  K3=S9J!lKNi 

%50  i  4i    TiQ2U)  Ti*         Flag4  M 

9660  i  42    Ti02:'s>  Ti4         Flaa4  M 

967(1  !  43    hN(i)  Tin         Tinflag  Khn 

9680  '  44    Ti'N(s)  Tin         Tinflaq  Ktin 

9690  i 

9700  ! 

9710  !  Lcxf lag  or  Lflaq  indica'es  presence  of  roft'nnoos  liquid  ohase 

o,120  !  AtO.  Ah,  At2.  At3.  At4  for  activity  hdeal  solutisn) 

«730  '  ?f  Set,  Ti0.  Ti203    h305,  T1Q2 

9?40  I  prr  :  activity  of  carbon 

975!  !  X  =  SQP CO?) 

9760  !  y  ^  S8R(H2) 

9770  i  Z  -  S8R(N2) 

9780  i 

'790  '      Variables  used  is  Newton's  netnod  calculation: 

9800  '  No  rtetai  or  o-ioe  X,Ti 

9B10  !  One  oxide  only/netal  only      <  (Ti  if  hN  *  i  solid  oxide,  no  wetal 

9820  '  Two  oxides/rtetal  +  1  oxide    - 

'8-i0  !  V  all  cases,  Z  and  Ace  are  *he  last  two  variables  wh*n  Gporonriatt 

984  0  END 
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KEY  0 

Check=(Chec)i=fi) 

-Execute 

KEY  i 

TRACE  VARIABLES  Yn(») 
-Execute 

KEY  2 

CONT  Ex 
-Execute 

KEY15 

-Clear  line 
SCRATCH 

KEY  4 

CONT  Temp 
-Execute 

KEY  5 

ftt$="Y" 

-Execute 

KEY  6 

TRACE  PAUSE  S 
-Execute 


KEY  7 


-Continue 


KEY  8 

Ht=.4*(Dlt=.i)  +  .i 
-Execute 

KEY  ? 

-Clear  line 
LOAD 

KEYifl 

-Clear  line 
SAVE 

KEYli 

-Clear  line 
STORE 

KEY12 

-Clear  line 
EDIT 

KEY13 

-Clear  line 
EDIT  LINE 

KEY14 

-Clear  line 
LIST 

KEY15 

-Clear  lme 
SCRATCH 

KEr23 
-11750 
-Continue 
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